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Abstract 

The computation of Feynman integrals in massive higher order perturbative calculations 
in renormalizable Quantum Field Theories requires extensions of multiply nested harmonic 
sums, which can be generated as real representations by Mellin transforms of Poincare- 
iterated integrals including denominators of higher cyclotomic polynomials. We derive the 
cyclotomic harmonic polylogarithms and harmonic sums and study their algebraic and 
structural relations. The analytic continuation of cyclotomic harmonic sums to complex 
values of N is performed using analytic representations. We also consider special values 
of the cyclotomic harmonic polylogarithms at argument x = 1, resp., for the cyclotomic 
harmonic sums at N — > 00, which are related to colored multiple zeta values, deriving 
various of their relations, based on the stuffle and shuffle algebras and three multiple ar- 
gument relations. We also consider infinite generalized nested harmonic sums at roots of 
unity which are related to the infinite cyclotomic harmonic sums. Basis representations are 
derived for weight w = 1,2 sums up to cyclotomy I = 20. 



Dedicated to Martinus Veltman on the occasion of his 80th birthday. 



1 Introduction 



The analytic calculation of Feynman integrals requires the complete understanding of the asso- 
ciated mathematical structures at a given loop level. By the pioneering work in Refs. [lj|2] this 
has been thoroughly achieved for one loop integrals occurring in renormalizable Quantum Field 
Theories. The corresponding complete framework in case of higher order calculations is, how- 
ever, not yet available. Which mathematical functions are of relevance there is revealed stepwise 
in specific higher order calculations. Here higher transcendental functions, like the generalized 
hypergeometric functions and their generalizations [3], play a central role. Their series expan- 
sion in the dimensional parameter [TJH] leads to nested infinite sums over products of digamma 
functions [5], cf. [5J. The nested harmonic sums [TUB] are defined by 

S b , a (N) = J2 {S,g lw ))k Sd{k) > S *( N ) = 1 > MieZ\{0}. (1.1) 

k=l 

They form a quasi-shuffle algebra [9J. Their values for N — > oo are the multiple zeta values 
resp. Euler-Zagier values [JU] defined by 

C M = lim S b 4N), b^l, 

iV— »oc 

see [IT]. In case of massless problems to 3-loop order the results for single scale quantities in 
Mellin space can be written by polynomial expressions in terms of S$(N) and (s with coefficients 
being from the rational function field Q(N), cf. e.g. [12]. In this context we consider the Mellin 
transform 

l 

N 



M[f(x)](N)= dxx N f(x). (1.2) 
Jo 

In most of the applications below we assume N G N + = N\{0}. 

In computations at even higher orders in the coupling constant and through finite mass effects, 
however, generalizations of the nested harmonic sums contribute, at least in intermediary results. 
One extension concerns the so-called generalized harmonic sums [13JH1] given by 

N C k 

S b , 3 (C,i;N) = ^2jsS^f]k), (1.3) 



k=l 



with b, di G N + ;£,£i G WL* = K\{0}. Known examples are related to the second index set 
& G {1, -1, 1/2, -1/2, 2, -2}, cf. lJ2J[15i[l6]. 

In case of single scale problems with two massive lines and mi = m 2 at 3-loop order summa- 
tion is also required over terms 



(1.4) 



(2k + l) n 



which is a special case of 



(1.5) 



(/ • k + m) n 



with l,m,n G N + . Note that we have deliberately chosen a real representation here, which is of 
practical importance in case of fast polynomial operations needed for solving nested summation 
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problems [T7] 0. Sums containing fractional terms m/l were considered in the context of colored 
harmonic sums in [20H22] . It is expected that sums of this kind and their iterations will occur in a 
wide class of massive calculations in higher order loop calculations in Quantum Electrodynamics, 
Quantum Chromodynamics, and other renormalizable Quantum Field Theories, in particular 
studying single distributions, but also for more variable differential distributions. Usually objects 
of this kind emerge first at intermediary steps and, at even higher orders, they occur in the final 
results. Therefore, these quantities have to be understood and methods have to be provided to 
perform these sums. This does not apply to integer values of the Mellin variable N only, but also 
to the analytic continuation of these quantities to iV £ C, since the experimental applications 
require the Mellin inversion into momentum-fraction space. 

We show that the single sums (11.51) and their nested iterations can be obtained from linear 
combinations of Mellin transforms of harmonic polylogarithms over an alphabet of letters con- 
taining x l /$k(x), < / < deg($fc(x)), where denotes the kth cyclotomic polynomial [23J. 
One may form words by Poincare-iterated integrals [25] over this alphabet, which leads to the 
cyclotomic harmonic polylogarithms fj, forming a shuffle algebra. This class extends the harmonic 
polylogarithms [26J. The Mellin transform of elements of has support x £ [0,1], but requires 
an extension in the Mellin variable N £ N, 

N^k-N, (1.6) 

where k denotes the index of This assumption allows to associate nested harmonic sums 

of the cyclotomic type by this Mellin transform. Special values are obtained by either the 
cyclotomic harmonic polylogarithms H^(x) at x = 1 or the associated nested harmonic sums for 
N — > oo. They extend the multiple zeta values and Euler-Zagier values, cf. Section [5j In the 
present paper we investigate relations and representations of these three classes of quantities. 
Special emphasis has been put on the class of nested sums f ll.4j) where all derived algorithms 
have been incorporated within the computer algebra package HarmonicSums [27] . 

The paper is organized as follows. In Section [2] we establish the connection between the cyclo- 
tomic harmonic sums and the cyclotomic harmonic polylogarithms through the Mellin transform, 
at modified argument kN . The basic properties of the cyclotomic harmonic polylogarithms are 
investigated in Section The cyclotomic harmonic polylogarithms obey a shuffle algebra. The 
nested sums at finite values of iV are studied in Section [U including their algebraic and structural 
relations, generalizing [28H31] . Here are also three multiple argument relations of interest. The 
analytic continuation of the new harmonic sums to complex values of N is presented based on 
recursion and asymptotic representations, similar to the case of the nested harmonic sum [2"9"H3"2"] . 
The representation of the cyclotomic harmonic sums requires to know their values at iV — > oo, 
which are equivalently given by the cyclotomic harmonic polylogarithms at x = 1. The set 
of special numbers spanning the Euler-Zagier and multiple zeta values [11] are extended. We 
study the cases of weight w = 1,2 up to cyclotomy I = 20 and derive the corresponding relations 
based on the stufffe and shuffle algebra and three multiple argument relations. Furthermore, 
we investigate the relations of the cyclotomic harmonic sums for iV — > oo extending for words 
resulting from the alphabet (±l) k /k, (±l) k /(2k + 1), cf. Section For the cyclotomic harmonic 
polylogarithms, harmonic sums and their values at N — > oo basis representations are derived. 
In Section [6] we study the relations of the infinite nested harmonic sums with numerators at lih 
root of unity, I < 20 for weight w = 1,2. Here, also the distribution relation, cf. [19], is considered 
beyond the relations mentioned before. We also consider a relation valid for finite generalized 

1 Complex representations are related to the so-called colored harmonic sums St,_s(p, f; N) = J2k=i ^^ct(f: 
with b,a % € N + ,p,n £ ufi 2 {exp[2ni(n/l)},n e {I,.-, I - 1}}, cf. (151113] . 
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harmonic sums with root numerator weights. Section [7| contains the conclusions. Some technical 
details are given in the Appendices. 



2 Basic Formalism 

We consider cyclotomic harmonic sums defined by 



N k 



S{ ai M,ci},...,{aiA,ci}( s ii — > s i'i N ) — ^ ^yr^5'{ a2ifeiC2 } ; ... ; { a;ibiiC! }(s2, s«; fci), So - 1, 

fcl=l ^ ^ 



(2.1) 



where a^, q G N+, 6, G N, Sj = ±1, cij > hi, the weight of this sum is defined by c\ + • — h q and 
{aj,6j,q} denote lists, not sets. One may generalize this case further allowing Sj G M*, [14J. Of 
special interest will be the infinite cyclotomic harmonic polylogarithms defined by 

0'{ai,6 X) ci},...,{o,,6,,q}(si,-..,Si) = Jim 5{ ai) 6 1 , Cl },...,{o,,6,,q}(si, iV) 

iV— >oo 

which diverge if C\ = 1, si = 1. Sometimes we will use the notation 

5'{ai,6i,c 1 },...,{a i ,6 i ,Ci}(' S l ; — > s h N) = S'{o 1 ,6i,a 1 ci},...,{oi,6j,s I ct}(-^') ( 2 - 2 ) 
cr {ai,bi,ci},...,{a i ,6;,c i }('Sl, S/) = 0"{ aii ;, ljSlCl } ) ... ) { aii b ijSiCi 

} ( 2 - 3 ) 

as shortcut below. 

For further considerations, we rely on the following procedure which transforms the sums to 
integrals. The denominators have the following integral representation 

(±1)k ,i 



(ak + b) c 
and the sum over k yields 



, dx x ak+b - 1 (±l) k (2.4) 

ak + b Jo v ' v ' 

;±D' /" p <fe r~ ^ r- ^ I **. l(±1) , i (2 . 5) 

JO x l x 2 Jo x e-l Jo 



^ (±1) V^-' = ^-' ^'_/ . (2.6) 



±a; a - 1 

k=l v 7 

This representation is applied to the innermost sum (a = ai,b = bi,c = q). One now may 
perform the next sum in the same way, provided a/_i|a/. If this is not the case, one transforms 
the integration variables in (12. 4[ 12. 5j) such that the next denominator can be generated, etc. In 
this way the sums ( 12. ip can be represented in terms of linear combinations of Poincare-iterated 
integrals. Evidently, the representation of the cyclotomic harmonic sum ( 12. ip in terms of a 
(properly regularized) Mellin transform will be related to the Mellin variable kN, with k the 
least common multiple of a\, ai. 

Let us illustrate the principle steps in case of the following example : 

N i k {-iy 

S{3,2,2 } , {2 , a ,i } (i, -i; N) = g £ (2IT1) • (2 - 7) 
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The first sum yields 

N 



Setting x = y 3 one obtains 



r 1 ? ,8 N (— v 6)k _ 1 



12/ ^ f^ifM"-' 



y & + 1 I 7 z J t 6 + 1 







y b + lj z J t b + l 



In general, the polynomials 



(1-7T) / f / ^^V^T^ (2-10) 



x a -l (2.11) 



in f)2.6p decompose in a product of cyclotomic polynomials, except for a = 1 for which the 
expression is $i(x), see Section 3. Moreover, the polynomials 

x a + l = - (2.12) 

x a - 1 

are either cyclotomic for a = 2™, n e N or decompose into products of cyclotomic polynomials 
in other cases, see Appendix A. All factors divide (x a ) 1 — 1, resp. (— x a ) 1 — 1. We remark that 
Eq. (I2.10p is not yet written in terms of a Mellin transform (II. 2p . To achieve this in an automatic 
fashion, the cyclotomic harmonic polylogarithms are introduced in Section 3. Furthermore, their 
special values at x — 1 contribute to which we turn in Section 5. 



3 Cyclotomic Harmonic Polylogarithms 



To account for the newly emerging sums (12.11) in perturbative calculations in Quantum Field 
Theory we introduce Poincare-iterated integrals over the alphabet 21 



21 



n r x l 

x) \ $ fc (x) 



k G N + ,0 < I < (p(k) } , 



(3.1) 



where $fc(a?) denotes the kth cyclotomic polynomial [23], and (p(k) denotes Euler's totient func- 
tion 



x n -l 



n ^ 

d\n,d<n 



(3.2) 



and the first cyclotomic polynomials are given by 



$1 ( x) 


= x — 1 




(3 3) 


$9 (x) 


= X + 1 




(3.4) 


<&i(x) 


= T 2 4- T 4- 1 






&d(x) 


= x 2 + 1 




(3.6) 


$e4x) 


= x A 4- x 3 + x 2 


+ X + 1 


(3 7) 


$fi(x) 


= x 2 — X + 1 




(3 8) 


x) 


= X 6 _|_ t; 5 _|_ x 4 


_i_ x 3 4- x 2 4- x 4- 1 


(3 9) 




- x 4 + 1 






$ 9 (x) 


= x 6 + x 3 + 1 




(3.11) 


$i (x) 


= x — X + X 


- X + 1 


(3.12) 


$11 (x) 


= x 10 + x 9 + x ! 


i + x 7 + x 6 + x 5 + x 4 + x 3 + x 2 + X + 1 


(3.13) 




= x 4 -x 2 + l, 


etc. 


(3.14) 



The alphabet 21 is an extension of the alphabet 

1 



2L 



1 

x' $i(x) ' $2 (^) 



(3.15) 



generating the usual harmonic polylogarithms [26J 0. As a shorthand notation we define the 
letters of 21 by 



i 

X 



x 



k E N+, Z G N, Z < . 



(3.16) 
(3.17) 



Here the labels k and Z form a double index, which always appear sat a common position. Mellin 
transformations associated to l/$4(x) and 1/§q(x) were discussed long ago in [36J, 



"(I 

with Stirling's /3-function [37] 

/3(x) 



dt 



t 



x-l 



t 2 + l 

1 



f3(x) 



dt t X 



x-l 



t-2 







X + 1 



t 2 - 1 + 1 



X 



(3.18) 
(3.19) 



(3.20) 



Integrals of this kind emerge also in particle physics problems in 



2 Note that we defined here the second letter by l/(x — 1) which differs in sign from the corresponding letter 
in . Numerical implementations were given in [33 i 34 . A few extensions of iterated integrals introduced in [215] 
based on linear denominator functions of different kind, which are used in quantum-field theoretic calculations, 
were made in [531155] . 
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We form the Poincare iterated integrals 

ct::!L( z ) = ^ $ (h, ■ ■ ■ , u = (o, . . . , o), (h, . . . , u = (o, . . . , o), 

Ctjz)= f dxf l C(x) iffc m ^0, 
Jo 

Ct±(z) = f dx C fc±>) if fa, ...,*») ^ (0, 0), 

Jo 

and C l s (z) denotes cyclotomic harmonic polylogarithms. They form a shuffle algebra [HIBJ] by 
multiplication 

C${z)-C%(z) = C%(z)inC%(z)= Yl C t^ ( 3 - 21 ) 

of M w elements at weight w, where M denotes the number of chosen letters from 21. The shuffle 
symbol LU implies all combinations of indices a^- € and fry G 6j leaving the order in both sets 
unchanged and the brackets [ ] pair the the upper and lower indices forming a unity, cf. (13.1 7J - 
The number of basis elements spanning the shuffle algebra are given by 

jV basic ( w ) = M d ,w>l (3.22) 

d\w 

according to the 1st Witt formula [H|H2]. Here \x denotes the Mobius function jl3]. The number 
of basic cyclotomic harmonic polylogarithms in dependence of w and M is given in Table 1. 



weight 


letters 


2 


3 


4 


5 


6 


7 


8 


1 


2 


3 


4 


5 


6 


7 


8 


2 


1 


3 


6 


10 


15 


21 


28 


3 


2 


8 


20 


40 


70 


112 


168 


4 


3 


18 


60 


150 


315 


588 


1008 


5 


6 


48 


204 


624 


1554 


3360 


6552 


6 


9 


116 


670 


2580 


7735 


19544 


43596 


7 


18 


312 


2340 


11160 


39990 


117648 


299592 


8 


30 


810 


8160 


48750 


209790 


729300 


2096640 



Table 1: Number of basic cyclotomic harmonic polylogarithms in dependence of the number of letters 
and weight. 



Now the cyclotomic harmonic sums (12. ip can be represented as Mellin transforms of cyclo- 
tomic harmonic polylogarithms. In the example (I2.10p the iterated integrals have to be rewritten. 
We express respective denominators in terms of products of cyclotomic polynomials, (IA.2HA.16]) . 



and perform partial fractioning in the respective integration variable, flA.27HA.38j) . In our con- 
crete example, the integrand of the y-integral in (12. 9 p has the representation 

■"' 1 [f°Ay)-fUy) + VUy)} ■ (3.23) 



y 6 + 1 3 
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With integration by parts one obtains the following Mellin transforms of argument 6N of cyclo- 
tomic harmonic polylogarithms C^'"'J^ n (x) weighted by the letters fj°(x) of the alphabet 21 : 



S. 



{3,2,2},{2,l,l}l i ' 



-1;N) = ±(4-ir)Jdxx 3 (x 6N -l)[6 + f°(x)-f°(x)-2f°(x) 

-ffa) -2f 6 (x) + ti(x)]C° (x) 

-2 f dxx* [(-l) N x 6N - 1] [3 - f!(x) - 2&{x) + 2f l2 {x)] C° (x) 
Jo 

~ K'4°(l) - Ol) + 2CSi(l)] f Q dxx* [(-Ifx™ - 1] 

x [3 - f!(x) - 2f° 12 (x) + 2fUx)] 

+1 ^ 1 rfxx 3 [(-1) V - 1] K; 4 °(x) - <f 2 (x) + 2C«%{x)} 

x [3-/?(x)-2tf a (x)+2y*(x)] . 

(3.24) 



The constants C^'Til) are discussed in Section 5. In particular one obtains 



a 



0,0 

0,4 



-C, 



(3.25) 



with C the Catalan number [33], and Cg' 12 and C*of 2 are linear combinations of ip'(l/12) and 
^'(5/12). The latter numbers reduce further, cf. Section 5. Note that in all cases in which neither 
fi(z) nor C^' fc 2 '"jfc ( z ) are present, the z-independent terms in [(±z k ) N — 1] can be integrated, 
since the other cyclotomic letters and C^|„"'^ (z) for fc > 1 are regular at z — 1. 

In general, linear combinations of the Mellin transforms 



M 



f&)-C b s (x) (IN) 



dxx lN (f?(x)r-C b d (x), uG{0,l}. 



(3.26) 



with I being the least common multiple of a 1; ...,a^, allow to represent all cyclotomic harmonic 
sums. Summarizing, we can write 



s {aiM , clh ..., {ak , h , Ck} (N) = J2 e J dx x lN (fZ:(x)) Un c%(x 



n=l 



(3.27) 



with e n G K and w n G {0, 1}; here e n is determined by polynomial expressions in terms of 
cyclotomic harmonic polylogarithms evaluated at 1 with rational coefficients. Within this trans- 
formation and in the following the letter 1/x plays a special role concerning the cyclotomic 
harmonic sums. We can exclude the case that the first letter in f!3.26j) is 1/x since it would just 
shift the Mellin index by one unit. In case of and related functions the +- regularization 



M 



/(*) 



x 



(iV) 



dx 



x 



N 



X 



(3.28) 



is applied to (I3.26p . Only in case that fi(x) = 1/ (x— 1) and C~(x) do not vanish in the limit x — > 1 
a +- function must occur. Iterations of cyclotomic letters f l k {x), k > 1,0 < I < k for k > 2 
have no singularities in x G [0, 1]. 
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We remark that the sketched transformation from cyclotomic harmonic sums to their Mellin 
transforms in terms of cyclotomic harmonic polylogarithms can be reversed, i.e., a given expres- 
sion in terms of Mellin transforms of cyclotomic harmonic polylogarithms can be expressed in 
terms of cyclotomic harmonic sums. The explicit algorithms have been worked out in details for 
the alphabet 

*' := {/ °, fl fl fl fl} = {-, ~^- y -±- y jfrlca, (3.29) 

1 J {X <Pl(x) <P2(£) ®4{X) <&4{X) J 

which allows one to express the cyclotomic harmonic sums 5{ai,bi,ci)v..,{aj,&i,cj}( s i) sf, N) with 
a» G {1,2}, bi G {0,1} and q G N + in terms of Mellin transforms of cyclotomic harmonic 
polylogarithms in both directions; the implementation is available within the HarmonicSums 
package [27]. This transformation will be used, e.g., in Section 4.1. 

The Mellin transforms (I3.26[ 13. 28j) obey difference equations of order I in N, which can be 
used to define these functions specifying respective initial values for I moments. In the following 
we illustrate this for the words ^'/^(x). One obtains 

c n ,kMN + n - = T7-T , (3.30) 



where $fc(x) is given by 



N 

n=0 



$ k {x) = Y,c n ,kX n . (3.31) 

n=0 



Here we define 



f 1 x N - 1 

<M0,iV) = / dx (3.32) 

Jo x ~ l 

MhN) = [ dxx N f k (x), <j> k (N)=<j> k (p,N), k>2 (3.33) 
Jo 

<f> k (l,N) + = [ dx[x N -l]f k (x), MN) + =M0,N) + . (3.34) 



For k < 105 for all coefficients c n>k G { — 1,0,1} holds [IS]. One derives the following first 
order difference equations f 1 3 . 3 j) for Mellin transforms associated to the lowest order cyclotomic 
polynomials : 

Ml,N+l)-M',N) = N+ ) + 1 (3.35) 

2 (/, N + 1) + M', N) = N+ ) + 1 (3.36) 

M l,N + 3 )- M ,,N) = - (JV + ; + 1) 1 (JV + ; + 2) (3.37) 

^ 4 ((,JV + 2)+^ 4 (i,JV) = N+ \ + 1 (3-38) 

M l,N + S)-M',N) = - (N + l + 1 ' {N + l + 2) (3.39) 

M l,N + 3)+M ,,N) = ^jtL + W— (3.40) 
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(AT + / + l)(JV + / + 2) v ' 

Ml,N + 4)-MKN) = N+ \ + 1 (3-42) 

Mt,N + 9)- M ',N) = - (7y + ; + 1) 3 (jv + ; + 4) (3-43) 

+ *0 = (iY+ 3 ^ i ;' ; ; + + 2 ; + 2) (3-44) 

M l,N + U)- M ,,N) = (jv + i + 1) 1 (Af + ; + 2) (3.45) 

+ = (Jf+ 3 /^ + + 2 l + 2) ' ^ (3 - 46) 

Together with the corresponding initial values, these recurrence relations enable one to compute 
efficiently the values for N. In particular, due the special form of the recurrences, we get explicit 
representations in terms of finite sums. E.g., for <J)q(1,N), we obtain 

M , a*) = (-!)» (± gJzM+fl^ _ + W , o)) , (3.47) 

and 



6 (Z,3iV + l) = 6 (/ + l,3AT) (3 
6 (/,3iV + 2 = 6 (Z + 2,3iV). (3.49) 



Looking at 



lim 6M, N) = lim / x^'— — dx = lim / a^+'fl + x - x 3 - x 4 + . . . )dx = 

f-Kx> JV-voo J X 2 — X — 1 N^oo J 

*.(,,«» = -£ ("DW^O-S] 



A' 

shows that 



^ (3z + Z - 2)(3i + Z — 1) ' 

Completely analogously, all the other functions (I3.35fl3"~4"t3j) can be written in such a sum repre- 
sentation where the constants are the infinite versions of it multiplied with a minus sign. Note 
that these constants can be written as a linear combination of the infinite sums 



0{a,b,s} = 



s k 



k=l 



ak + b 



with s G {—1, 1} and a, b G N with a / 0. The relevant values for this article will be worked out 
explicitly in the Section 5. 

The cpk(l,N) + functions are preferred if the recurrences (13.35143^61) form telescoping equa- 
tions. In this case, in particular if k is odd, the riV) for properly chosen r can be related 
to sums without any extra constant, e.g., 



N 



5 (O5iV) + = -|: (5 . + t _ 4)(5 . + t _ 3) . (3.50) 
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i>i and 02 can be related to the single cyclotomic harmonic sums at weight w 
1 (/,iV) + = S 1 (N + l), 
fa(l,N) = (-If '\S-!(N + I) -ln(2) + S-!(l) . 



1 as follows 

(3.51) 
(3.52) 



For later considerations we use S±(N) and S-i(N) instead of </>i(iV) and <f>2(N). In the following 
Section the single cyclotomic harmonic sums of weight w = 1 are expressed in terms of the Mellin 
transforms 4>k{l, r N) for properly chosen r. In addition, the cyclotomic harmonic sums of higher 
weight and depth will be discussed. 



4 Cyclotomic Harmonic Sums 

We consider the extension of the finite nested harmonic sums (TIE] to those generated by the 
cyclotomic harmonic polylogarithms discussed in Section 3. First the single cyclotomic harmonic 
sums are considered and explicit representations are given. We derive their analytic continuation 
to complex values of N. Next the algebraic, differential and three multiple argument relations of 
the cyclotomic harmonic sums are discussed. These relations are used to represent these sums 
over suitable bases. Finally we consider the nested sums over the alphabet {(±l) fc / k, (±l) fc / (2k+ 
1)} to higher weight deriving explicit relations up to w = 5. 



4.1 The Single Sums 

The single cyclotomic harmonic sums are given by 

^2 (±1)k 

k=0 



(l-k + to) 7 



(4.1) 



Here N is either an even or an odd integer. In case one needs representations only for iV G N 
the following representations hold : 



N (_-i\k 
^— ' l ■ k + to 

gtj (~l) fc 

' I • k + to 

fc=0 



N 



E 

.fe=0 
' N 

E 

,fc=0 



1 



(2/) ■ k + to (2/) • k + to + I 
1 1 



(2/) • k + to (2/) -k + m + l 



(21) N + l + m 



(4.2) 



(4.3) 



with iV = 2N. However, one is interested in relations for general values of N, since for nested 
sums more and more cases have to be distinguished. The single sums can be expressed in terms 
of the Mellin transforms 4>k{l, kN), (I3.32f[3~33l . Up to w = 6 one obtains : 



- 1 

^l + 2k 

k=i 

f (-D- 

^ l + 2k 

k=i 

- 1 

T — 
^ i + 



2N Si(N) 



2N+1 2 
1 



Si (2 iV), 



2N+1 



4 (2iV) 



+ C {2,1,-1}, 



k=l 



3k 



3N+1 6 2 u ' 2 Viy ; + ' 



(4.4) 
(4.5) 
(4.6) 
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N 

E 

k=i 



l + 3k 



2 + 3k 



N 1 

fe=i 

y (-i)' 

f--2 + 3t 

fe=l 



Er^ 



1 + 4A; 

(~l) fc 
1 + 4A; 

E — 

^3 + 4k 



k=i 

N 

E 

fc=i 



fe=i 



3 + 4k 



E 

fe=i 



E — 



k=i 



N 

E 

fc=i 



1 + 5/c 



(~l) fc 
1 + 5A; 



" 1 

9~r 



fe=i 



E 

fc=i 



(~l) fc 
2 + 5k 



N 1 

Q i 



k = l 



N 

E 

k=i 



(~l) fc 
3 + 5A; 



^(TV) - ^(STV) + (-1)* 



37V 



1 



37V + 1 2 



o0e(3iV) 



^{1,0-1} + (7 {3 ,1-1}, 



(4.7) 
(4.8) 



^-i(iV) + ^(STV) + (-1)" 



+ g<7{l,0,-l} + <7{3,l,-l} + 2' 



1 



3iV + 2 2 



o0e(3iV) 



27V 1 



= (-1) 



— - lSl (2N) + -*(4JV) - -0 4 (47V) + -. {2!l ,_ 1} + 
1 



A? 



47V + 1 



*(4iV) 



+ <7{4,1,-1}, 



(4.9) 
, (4.10) 

(4.11) 



107V 



3(4JV + 3) - i S ' (2jV) + ^ + 2^ (47V) - r<^-» - 2(57+8)' 



=(-1) 



A? 



47V + 3 



8 (2,47V) 



+ ^{4,3,-1}, 



(4.12) 
(4.13) 



- ^ + ^(57V) - ^ 5 (57V) + - ^> 5 (1,57V) + - ^ 5 (2,57V) + , 



(4.14) 



:-i) 



IV 



-\<t> w {$N) + I^ 10 (l, 57V) _ l0 lo(2 , 57V) + _JL_ 



+ ^S-i(JV) - ^-i(57V) + ^ {1A _i } + <7 {5l i,-i}, 



(4.15) 



5iV 



Si (TV) , 1 



2 + 5A; 2(57V + 2) 20 



+ -S^N) + -0 5 (5iV) + - -0 5 (l,5iV) + - -0 5 (2,57V)+, 



-1) 



N 



-\<h>(5N) ~ ^io(l, 57V) + ^io(2, 57V) + 



(4.16) 
(4.17) 



^S^(N) + \s^(5N) - ^ {1 , ,_ 1} + (7 {B>2> _ 1} , 



57V 



Si (TV) , 1 



1 



1 



3 + 5k 3(57V + 3) 20 



+ -5i(57V) + -0 5 (57V) + + -0 5 (1,57V) + - -0 5 (2,57V) 



+ ! 



(4.18) 



=(-1) 



N 



i0 lo (57V) - ±0 1O (1, 57V) - lfc (2, 57V) + 



+ ^ S ~l( N ) ~ ^-l(57V) + ^{1,0,-1} + ^{5,3,-1}, 

12 



(4.19) 



N I 

^ A 4- 



fc=l 



E 

fc=i 



4 + 5A; 



(~l) fc 
4 + 5A; 



* 1 



fc=l 



jv , 

E{ 



+ 6A; 



-1)' 



+ 6A; 

E — 



fc=i 

N 



k=l 



N 

E 

fc=i 



5 + 6/c 



(~l) fc 
5 + 6A; 



+ ^5(1, 5JV)+ + |0 5 (2, 5JV)+, 



(4.20) 



-ifc„(5JV) + ~0 1O (1, 5N) - ^ 10 (2, 5iV) + 



113 7 

- ^Tj 5 - 1 ^) + 4 5 -l( 5iV ) + 5 a {lA-l} + 0-{5,l -1} - ^{5,2,-1} + 0-{5,3 -1} + ^ 



(4.21) 



()A + - 7 5 i( 2Ar ) - iS^N) + ^(6iV) - ^0 3 (3iV) + - ^0 3 (6iV) + , 



6N + 1 



12 



6 



-0i 2 (6iV) + i0 4 (2iV) + 



+ C{6,1 -1}, 



(4.22) 
(4.23) 



6N 



+ ^i(A0 _ ^ _ l s , 3N ) + I 5 ( 6A n 
5(6iV + 5) 12 6 V 7 4 v 7 2 v ; 

+ ^0 3 (3iV) + + ^ 3 (6iV) + , 



(4.24) 



-1 



0i 2 (6iV) - 2 -Mm - MQN) + — ±— 

3 bJS + 5 



+ ^{2,1-1} ~ <7{<&,i,-i} + j^- 



(4.25) 

Taking the sum representations such as (13.501) and (13.471) given by the recurrence relations H3.35H 
13.461) and the corresponding initial values </>&(/, 0) (which are expressible by the infinite sums 
C{a,6,±i}) cf. Section 5), we used the summation package Sigma [UJ to perform this transforma- 
tion. As a consequence, the single cyclotomic harmonic sums can be expressed in terms of the 
sums 



S-i(iV), 5_ x (3iV), S^(5N), S 1 (N), S 1 (2N), S 1 {W), S 1 (AN), S 1 {5N),S 1 {6N), 

3 (3iV)+, 3 (6iV) + , 4 (2iV), 4 (4iV), 4 (6iV), 5 (5iV) + , 5 (1, 5JV)+, 

5 (2, 5A0+, 6 (3iV), 8 (4iV), 8 (2, 4JV), 0io(5iV), 1O (1, 5JV), 0io(2, 5iV), <M6iV). 

In particular, in the way how this construction is carried out it follows by the summation theory 
of [17] that the sequences produced by these sums form an algebraic independent basis over 
the ring of sequences generated by the elements from R(iV)[(— 1) N ]. Note that we exploited the 
algebraic relations (|5.118ff5\129p which implies that only the following constants 



C{1,0,-1}> cr {2,l,-l} 5 °"{3,1 ,-l}> C {4,1,-1}; C{4,3,-1}> C{5,1,-1}, °"{5,2 ,-1}, 0"{5,3-l}°"{6,l,-l} 

from R appear in the representation found for the single cyclotomic harmonic sums with weight 
w = 1. 

Due to (I3.33[ I3.34j) the functions <3?fc(/, N), k > 1 can be represented as factorial series [36JH6]. 
Therefore they are meromorphic functions in iV with poles at —n, jigN. This also applies to 
<f>k(Q,N), (I3.32p . The latter function grows oc ln(iV) for N — > oo, |arg(iV)| < ir. The recursion 
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relations (I3.30p allow one to shift 4>k(l,N) in N — > N + 1. To represent a function (f>k(l,N) for 
JVeC one needs to know its asymptotic representation in addition. 

It is given in analytic form in terms of series involving the Stirling numbers of the 2nd 
kind [36j[37]. The corresponding representations read : 

<M0,iV) ~ 7 + ln(iV) + ^-^ + -^--i-+ ' 



2N UN 2 120iV 4 252iV 6 2A0N 8 

691 0(±) (4,6) 



132iV 10 32760iV 12 V^ 13 , 

_1 1_ _L 1_ 17 _ 31 691 / 1 \ 

02(, ) ~ 4iV2 + 8iV4 4^ + 167V 8 iiV 1 ^ + 8^ + V^ 13 / 

, , rt *rN 1 2 2 14 1618 3694 1 \ /..ox 

~ W-^ + 3^-^ + 2W-3Am + °(iv5j (428) 
, . _ 1 1 5 61 1385 50521 A / 1 \ . , 



M0,N) 



1 1 1 1 31 67 109 361 3779 

^ ' ' ~ 5iV + 5iV2~ 5iV3~ iV4 + 25iV5 + 5iV6 ~ 5iv^ ~ iV 8 + 5^ 

412751 214093 1150921 n / 1 \ ,„ , 

+ 25^"liVir-^V^ +0 (ivHj (4 ' 30) 

, i 1 2 22 602 30742 2523002 n / 1 \ . , 

, , n ArN 1 2 16 12 56 3900 20296 5796 

67(0 N) ~ 1 1 1 

r7K ' ; 7AT 7 A^ 2 7N 4 7iV 5 iV 6 49N 7 7N 8 AT 9 

1809992 4582500 35282968 n / 1 \ 

+^ir + ^n ^ + o(^J (4.32) 

. , 1 1 1 3 11 57 361 2763 24611 250737 

6a(0,N) ~ 1 1 1 1 

V8K ' ; 2iV 2N 2 2iV 3 2N 4 2iV 5 2N 6 2N 7 2N 8 2N 9 

2873041 36581523 512343611 / 1 \ 

+^NTo 2W~i 2N^ + \N^J (4 ' 33) 

1 2 2 28 34 1172 1862 101428 207394 

09(0, ) ~ 3N + 3W~ 3m~3N I + 3N E + ~3W ~~ 3iV^ ~ 3iV 8 + 3N 9 

14999012 37996022 3386034628 / 1 \ 

+ 3iV 10 SiV^ 3A^ + \N^J ( ' ^ 

1 1 5 17 151 871 11465 92777 1626151 

0io(O, ) ~ N + Jp~Jp~m + W + W~H ; P iV 8 " + N 9 

16922791 370714025 4715323337 1 \ , 



1 4 6 56 282 3064 26646 382616 



UN 11N 2 11N 3 UN 4 11N 5 llN e 11N 7 11N 8 
4592442 7618184 13945859346 28200213176 n / 1 \ 

+ ,nn + ini »m ,»m + ° Urn ( 4 ' 36 ) 



hat 9 atio miv 11 HA 7 " 12 Va^ 13 



1 1 7 23 305 1681 33367 257543 6815585 
™ ' ^ ~ N + N 2 ~~Jp~m + N^ + ~~W W N*~ + N 9 

67637281 2237423527 27138236663 „ / 1 \ tl n . 
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, . , T , 1 2 86 338 12094 690866 / 1 , 

_ . . .280 25866752 ^ / 1 , 

^s(2,iV) ~ - - + _-_ + — __ + of — ] (1.39) 



1 2 


40 

+ iV5 


1952 


1 

+ - 


2iV ~ iV3 


N 7 


1 1 

iV + N 2 


5 

N 3 


17 

AT4 + 


151 



.... . . 871 11465 92777 1626151 

16922791 370714025 4715323337 1 \ , 

+~ JfS jvH— +0(jplj M>) 

, . Art 1 6 186 14166 2009946 458225526 A / 1 \ 
M2,iV) ~ ___ + ___ + _ _ + (_j . (4.41) 

The numerical accuracy of the asymptotic representations at a given suitably large value of TV 
lowers with growing k, i.e., one has to choose larger values of N correspondingly to apply the 
asymptotic formulae. The above expansions can easily be extended to higher inverse powers of 
N. The recursion for <f>k(l, kN) is given in fl3.30p . resp. f)3.351 - [3~4^]) in a more compact form. 
Due to these for any iV e C at which (j)^{l,N) is analytic one may map (j>k(l,N) to values 
\N\ ^> 1, arg(iV) < 7r and use the asymptotic representations. 

The single cyclotomic harmonic sums of higher weight obey the representations 

V——— = ) } - / dx In"" 1 (a;) x m ~ l % , , m < I (4.42) 
f^ Q (lk + m) n (n-l)\J KJ x l -l K ' 

^z} (—i) k (—~\) n r 1 (—r l ) N — '\ 

Y } , = / > / dx \n n -\x) x m ~ l { - , m<l, (4.43) 

^ Q {lk + m) n (n-l)\J KJ x l + l v ' 

for l,m,n G N+, which may be expressed in terms of cyclotomic letters again. We note that 

(> " 1 dxx m - l f{x)=f dxx™- 1 ln n (x) f(x) . (4.44) 



o 



Therefore, (14.421 14~4"3)) can be expressed by the corresponding derivatives of </>&(n, IN) for iV and 
corresponding constants. 

4.2 Cyclotomic harmonic polylogarithms at x = 1 

We consider the cyclotomic harmonic polylogarithm C c 'i(x). Its value at x = 1 is given by 

Cfjl) = f dxf*(x)C<j(x) . (4.45) 



a,b 

J 







Let = x /$fe(x), with / < deg($^(x)) and n be the smallest integer such that (x) | (x n — 1), 
with 

iiy ' (4.46) 



$ fc (x) x n - 1 
Since x G [0, 1], the representation 



x 



-x l H$k(x)J2 xjn = J2 a ™^ xm I2 xJn > a m(l)eZ (4.47) 



KV 1 k j=0 m=0 j=0 
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holds. Thus we get 



f 1 x l w °° 1 

u v ' m=0 j=0 

representing the value of the depth d=l cyclotomic harmonic polylogarithms at linear 
combination of the corresponding infinite cyclotomic sums, see Section 
We return now to Eq. (14.451) . Integration by parts yields 

00 w pi jn+m+l 

C$1) = C a c (l)cf(l) -£$>™ / dx . X C({ X ) . (4.49) 

Since the Mellin transform of a cyclotomic harmonic polylogarithm can be represented by a 
linear combination of finite cyclotomic harmonic sums, the weighted infinite sum of the latter 
ones in (I4.49p is thus given as a polynomial of infinite cyclotomic harmonic sums, Section 5. 



4.3 Relations of Cyclotomic Harmonic Sums 



The cyclotomic harmonic sums obey differentiation relations, cf. [7 
their quasi-shuffie algebra, cf. [9], and multiple argument relations 



stuffle relations due to 



4.3.1 Differentiation 



As has been illustrated in Section 3, cyclotomic harmonic sums can be represented as linear 
combinations of Mellin transforms f!3.27j) of cyclotomic harmonic polylogarithms. Based on this 
representation, the differentiation of these sums is defined by 

^S {ai , 6l , cl} ,..., {afciW (iV) = J> n / dx x lN l m ln m (x)(/f : (x))-C^(x). (4.50) 

n=l ^° 

The product \n m (x)C$? (x) = rn!CQ'"*'gC^! l (a;) with m consecutive zeros may be transformed into 
a linear combination of cyclotomic harmonic polylogarithms using the shuffle relation (13.211) . 
Finally, using the inverse Mellin transform, the derivative (14.501) of a cyclotomic harmonic sum 
w.r.t. N is given as a polynomial expression in terms of cyclotomic harmonic sums and cyclotomic 
harmonic polylogarithms at x = 1. Together with the previous section, the derivative (I4.50P can 
be expressed as a polynomial expression with rational coefficients in terms of cyclotomic harmonic 
sums and their values at iV — > 00. The corresponding relations are denoted by (D). Further 
details on cyclotomic harmonic sums at infinity are given in Section 5. 

A given finite cyclotomic harmonic sum is determined for iV 6 C by its asymptotic represen- 
tation and the corresponding recursion from iV — > (N — 1). Both the asymptotic representation 
and the recursion can be easily differentiated analytically. Therfore any differentiation of a cy- 
clotomic harmonic sum w.r.t. iV is closely related to the original sum. For this reason one may 
collect these derivatives in classes 

S[Zuc l} ,..M kM (N) = J^ { „ lM} ,..., Wrt} (iV);n G n| . (4.51) 

3 For the relations given in this Section we mostly present the results, giving for a few cases the proofs in 
Appendix B. The other proofs proceed in a similar manner. 
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4.3.2 Stuffle Algebra 

To derive the stuffle relations, cf. [HJ|28], we consider the product of two denominator terms. 
For ai, a 2 , bi, 62, ci, c 2 , i G N they are given by 

C1 + C2— J — 1\ 1 



+ ( mf( ivY Ci + C2 ~- 7 '~ 1> \ a ^" Ja ^ 1 



and 



(ai2 + 6i) c i(a 2 « + & 2 ) C2 \a 2 / (aii + bi) c ^ 



> c 2 

ai 



(4.52) 



(4.53) 



if ai& 2 = a 2 6i. 

The product of two cyclotomic harmonic sums has the following representation. Let 
d, bi, di, ej, k,l,n G N + and Cj, /; GZ \ {0}. If aiei 7^ di&i one has 

5'{a 1 ,bi,ci},...,{a fc ,6 fc , Cft }(^)5'{ (ilieii/l } i ... i { diie! j,}(n) = 

n / si 

S!g n (Cl) „ , .x „ / .x 

/ , _|_ ^ 1 ^|ci| { a 2^2,c 2 },.--,{afc,''fc,Cfc}l Z / D {di,ei,/i},...,{d i ,ei,/ ; } v^J 

^ sign (fx) 1 

+ M i _|_ e )|/ 1 |' b {»i.6i.ci}.-">{afc,bfc,c fc }W<3{d2,e2,/2},...,{<i,,e i ,/ l }l«J 



- v f-iv-i V(-ivY |cxl + 1/11 " J " M 1 

) f^; } \ l/il-l / aiei - di&x (ai* + 



i=i 

X 'S'{a2,fe2,C2},.--,{ a fe: b fc: C fe}(^)'^'{ rf 2, e2,/2},..,M,e i ,/ i }(«), (4-54) 



resp. for a±ei = d\b\ one has 

5'{a 1 ,6 1 , Cl },...,{a fc ,6 fc ,c fc }( n )'S'{di,ei,/i},...,{ a ! i , e! ,/ i }(^) 



^ sign(ci) 1 

2^ ( a i + b ^|c 1 | ^{ a 2' fe 2- c 2}v--.K^^fc}l Z PK,eiJi},...,{rf i , e; ,/ i }l^ 

^ sign (A)' M 

+ 1^ (dtf + ei ) |yi I ti {^M,ci},...,{a k ,b k ,c k }{ i )^{d2,e2j2},...,{di,e l J l }^) 

<h\ 1/11 ^ ( sign sign 



' i=i 

^ S {a2,b 2 ,c 2 },...,{a k ,b k ,c k }{i)S{d 2: e 2 J 2 },...,{d l ,e u f l }{i) • (4.55) 



Subsequently, the relations given by (I4.54[ [435]) are denoted by (A). 
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4.3.3 Synchronization 

A first multiple argument relation is implied as follows. Let a, b, k G N, c G Z\ {0}, k > 2. Then 



fc-i 



5 {aAc} (fc ■ iV) = J2 ^gn (c) 4 5 {fc . ajfe _ a . JjSign(c)ft|c|} (TV) . (4.56) 

i=0 

For aj, 6j, m, fc G N, q G Z \ {0}, k > 2, the general cyclotomic harmonic sums obey 

^{a m ,bm,Cm},{am-l,bm-l,C m -l},— ,{oi,6i,ci} ■ -A) — 

5{a m _i,6 m _i,c m _i},...,{ai ) 6i ) ci} (fc ' j ~ i) sign (c m ) fc ' J ~* 



EE 



(4.57) 



Repeated application of (14.561 I4.57j) allows to represent cyclotomic harmonic sums of argument 
kn by those of argument n. Subsequently, the resulting relations are denoted by (M). 

4.3.4 Duplication Relations 

The usual duplication relation [IT] holds also for cyclotomic harmonic sums : 

E s < 

},...,{ai,6i,±ci} (2A) — 2 <S'{2o m ,6 m ,cm},...,{2o ll 6i,ci 

} (A) . (4.58) 

± 

The proof of (I4.58P is given in Appendix B. 

Similar to (14.581) one obtains a second duplication relation : 

• diS{ 

} (2-A) 

die{-i,i} 

— 2 5'{2o m ,6 m — a m ,c TO },...,{2oi,6i— oi,ci} (A) , 

(4.59) 

its proof is given in Appendix B. The resulting algebraic relations of (I4.58P and (14.591) are denoted 
by (Hi) and (H2), respectively. We remark that more general relations arise for generalized 
cyclotomic harmonic sums presented in Section 6, (16.121 IBTTBl . 

4.4 Sums of Higher Depth and Weight 

As an example we consider the class of cyclotomic sums, which occur in the physical application 
mentioned in Section 1. They are given by iteration of the summands 

1 (-l) fe 1 (-l) k 

— — — — (4 60) 

jfe'i' fck ' (2k + 1) 1 ^ (2k + 1) 1 * y j 

with N > k > 1, k,k £ N + . As pointed out earlier, this class of cyclotomic harmonic sums can 
be expressed by Mellin transforms in terms of cyclotomic harmonic polylogarithms generated by 
the alphabet fl3T29l . 

We apply the relations in Section 4.2 to derive the corresponding bases for given weight w. 
The cyclotomic harmonic sums can be represented over the corresponding bases. Applying the 
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relations using computer algebra we find the pattern for the number of basis elements up to w 
= 5 given in Table 2. 



w 


N s 


Hi 


Hi, H 2 


Hi,M 


H U H 2 ,M 


D 


Hi,H 2 ,M,D 


A 


H U H 2 ,M,A 


A, D 


all 


1 


4 


3 


3 


2 


2 


4 


2 


4 


2 


4 


2 


2 


20 


18 


17 


16 


15 


16 


13 


10 


8 


6 


6 


3 


100 


96 


93 


92 


89 


80 


74 


40 


35 


30 


27 


4 


500 


492 


485 


484 


477 


400 


388 


150 


142 


110 


107 


5 


2500 


2484 


2469 


2468 


2453 


2000 


1976 


624 


607 


474 


465 



Table 2: Reduction of the number of cyclotomic harmonic sums Ns over the elements (I4.60P at given weight 
w by applying the three multiple argument relations (H\, H2, M), differentiation w.r.t. to the external sum 
index N, (D), and the algebraic relations (A). A sequence of symbols corresponds to the combination of 
these relations. 



Due to the arguments given in Section 4.3.1 above, for any occurrence of a differential operator 
(d m /dN m )S ai fi 1 , Cl , a>k, bk, Ck(N) only one representative is counted. 
The total number of sums of weight w for the alphabet f)4.60p is 

N s (w) = 4 ■ . (4.61) 

We mention that the single application of any of the multiple argument relations (M, Hi, H 2 ) 
leads to the same number of basis sums 

N Hl = N H2 = N M . (4.62) 

This also applies to the combinations H\M and H 2 M 

N Hi ,m = iWr- (4.63) 

Explicit counting relations for the number of basis elements given Table 2 can be derived : 

N A (w) = ^EKl) 5 " (4 ' 64) 

d\w 

N D (w) = N s (w) - N s (w - 1) = 16 • 5 W ~ 2 (4.65) 

N Hl (w) = N s (w) - 2 W - 1 = 4 • 5" 1 - 1 - T"- 1 (4.66) 

N HlH2 (w) = N s (w) - (2 ■ 2™- 1 - 1) = 4- 5™- 1 - (2- 2^-1) (4.67) 

N HlM (w) = N s (w) - 2 ■ T~ x = 4 ■ - 2 ■ 2 W ^ (4.68) 

N Hi h 2 mH = N s (w)-(3-2 w - i-l)=4-5 w - 1 -(3-2— 1 - 1) (4.69) 
N AD (w) = N A {w) - N A (w - 1) 



d\w d\w-l v 7 

d\w \ d\w I 

N DHiH2 mM = N HlH2M (w)-N HlH2M (w-l) = lQ-5 w - 2 -3-2 w - 2 (4.72) 
N ADHi h 2 m{w) = N AHlH2M (w) - N AHlH , 2M (w - 1) (4.73) 



d\w d\w-l v 7 



(5 d - 3 • 2 d ) . 

(4.74) 
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Here /i denotes the Mobius function 

The analytic continuation of the cyclotomic harmonic sums can be performed as outlined 
for the case of the single sums in Section 4.1. Their representation as a Mellin transform of 
the cyclotomic harmonic polylogarithms (I3.26P relates them to factorial series, except the case 
c\ — 1, si = 1, cf. ( 12. ip . If a sequential set of first indices q = 1, Sj = 1 occurs, one may reduce 
the corresponding sum algebraically to convergent sums, separating factors, cf. [5]. I.e., also 
in the general case the poles of the cyclotomic harmonic sums are located at — k, k e N. The 
recursion relations of the cyclotomic harmonic sums (12. ip imply the shift relations iV — > N + 1 
in a hierarchic manner, referring to the sums of lower depth. 

To accomplish the analytic continuation in N, the asymptotic representations of the cyclo- 
tomic harmonic sums have to be computed. Since the cyclotomic harmonic sums are represented 
over respective bases, only the asymptotic representations for the basis elements have to be de- 
rived. One way consists in using iterated integration by parts 



dx x N fix) = -j^t- - — / dx x N \xf'(x)] , (4.75) 

with f(x) a linear combination of cyclotomic polylogarithms. Here f(x) is conveniently expressed 
in a power series to which we turn now. 

We illustrate the principle steps considering the alphabet (I3.29P related to (I4.60p . cf. (13.161 
I3.17p . In general the cyclotomic harmonic polylogarithms Cg(x) over the alphabet ( I3.29P do not 
have a regular Taylor series expansion, cf. [26]. This is due to the effect that trailing zeroes in 
the index set may cause powers of ln(x). Hence the proper expansion is one in terms of both x 
and ln(x). For depth one and < x < 1 one obtains 

C °(x) = ln(x) (4.76) 



CS(x) = - E ^L = - E ^ + Eii L_ (4.77) 

1=1 8=1 1=1 

<*<*) - -Ey-E^ + E 1 ^ < 4 - 78 > 

i=l i=l i=l 

°° (_iy T 2i-i 



Cl(x) = (4.80) 

i=l " 

Let C\[x) be a cyclotomic harmonic polylogarithm with depth d. Assume that its power series 
expansion is of the form 

j=l i=l v 

with the index sets a and 6 according to the iteration of the letters (l4.76fT4T80|) and fij a corre- 
sponding index set of the cyclotomic harmonic sum, x e (0, 1), w £ N and Cj G Z. 

Then the expansion of the cyclotomic harmonic polylogarithms of depth d + 1 is obtained by 
using 



" x a'x 2 ^ 
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'{2,Cj ,cra},rij 



°1% X ) = EEf9u,.^n 5 {2, Cj ,-},% (*) + E E Wi. T7^*\ s t 

(4.83) 



^ ^ (2i + c, + 1) ~^>™^ w ^ ^ (2z + c ■ + 2) 

3=1 4=1 v J j = i j=i V J / 



3 = 1 4=1 v J ) j = i j=i V J / 

(4.84) 

W OO y 2i+ Cj + l 

^» = EE a + c. + n W^W (4-85) 

3 = 1 4=1 3 1 

w °° f_iy T 2i+Cj+2 

^S(^) = EE 2i + , + 2) V,-^(i (4-86) 

3 = 1 4=1 ^ 3 1 

Sample proofs of these relations are given in Appendix B. 

The analytic continuation of the cyclotomic harmonic sums at larger depths to iV G C is 
performed analogously to the case discussed in Section 4.1. Shifts parallel to the real axis are 
performed with the recurrence relation induced by (12.11) . The asymptotic relations for N — > 
oo, | arg(iV) | < 7r can be derived analytically to arbitrary precision. Examples are : 

, a ... ... ... /vr 2 1 



%l,2},{l,0-2} (N) ~4CfOjo(l) + 4C / o(l)C / o J o J o(l) - 4C / o J O i/ O i/ o(l) + ( — - — 
+ 



2 N 

I II 3 127 5 221 7 



N 2 12N 3 AN A 240iV 5 16iV 6 1344iV 7 64iV 8 

9 W(i) + (-W-U-^ + 27 



384(W 9 256iV 10 / f *>f° K ' 1 ; \IW 4 AN 5 64iV 6 

1 269 11 8699 \ „/ 1 \ , „ N 

+ +Ohm (4-87) 



32iV 7 256iV 8 32iV 9 102AN 10 J \N n y 

SW,2},{2,l,-2} (iV) ~ - y + CfoMCfcfafoQ) - C f o )f o )fhf o(l) + -L - ^ + 

3 127 5 221 7 367 

~ 16iV 4 960iV 5 ~ 64iV 6 + 5376iV 7 ~ 256iV 8 + 15360iV 9 

9 + ( _W i _ s + 25 « 



1024iV 10 v y l64iV 4 64iV 5 128iV 6 256iV 7 



77 221 1545 

+ 7777777 + 



10 



( 7f2 



1024iV 8 512iV 9 1024iV 
1 1 11 3 127 5 221 

+ 777 - 7777 + 777777 ~ 777777 + 77777777 ~ 777777 + 



4AT 4AT2 48AT3 16AT4 qqon$ 64iV 6 5376iV 7 
7 367 9 \ / 1 \ 

"2567V 8 + 1536(W 9 ~ 1024/V 10 J + V^V ' 

Given a cyclotomic harmonic sum with the iterative denominators (14.601) and given the number 
of desired terms, the corresponding expansion can be computed on demand by the HarmonicSums 
package. 
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5 Special Values 



The values of the cyclotomic harmonic polylogarithms at argument x — 1 and, related to it, 
the associated cyclotomic harmonic sums at N — > oo occur in various relations of the finite 
cyclotomic harmonic sums and the Mellin transforms of cyclotomic harmonic polylogarithms. 
In this Section we investigate their relations and basis representations. The infinite cyclotomic 
harmonic sums extend the Euler-Zagier and multiple zeta values [TT] and are related at lower 
weight and depth to other known special numbers. We first consider the single non-alternating 
and alternating sums up to cyclotomy I = 6 at general weight w. Next the relations of the infinite 
cyclotomic harmonic sums associated to the summands (I4.60p up to weight w = 6 are worked 
out. Finally we investigate the sums of weight w = 1 and 2 up to cyclotomy I = 20. 



5.1 Single Infinite Sums 

We consider the single sums of the type 

f (**)* (5 1) 

^ (Ik + m) n ' V ' } 

with l,m,n G N + ,l > m,n > 1. These sums are linearly related to the colored harmonic sums 
by 

OO I OO - / O ' \ 

£^X>rgp^p (5.2) 

and similar linear relations for the nested sums. We will address the latter sums in Section |6j 



5.1.1 Non-alternating Single Sums at w = 1 

For the non-alternating sums one obtains 



E 

fc=0 



Ik + m 



1 

1 L 



^0 - IE - lj) ( y 



m 



with 



OO j 

5^ k 

k=l 



(5.3) 



(5.4) 



denoting the divergence in (15. 3p . The corresponding sums, except cr can always be regularized. 

The reflection symmetry of the -^-function interchanging the arguments x and (1 — x), [36] . 
implies 



E 



fc=0 



1 



1 



7T , 




= — cot 


K n*) 


n 





nk + I nk + (n — /) 
The digamma-function at positive rational arguments obeys [551117] 

[(?-l)/2] 



(5.5) 



-Je - hi(2g) - - cot ( — 



= -n(7 S + ln(n))-^V Qj • 



/ 2nkp\ 



q 



) 



In 



sm 



q 



(5.6) 
(5.7) 
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Eq. ( 15. 6 p is used to remove dependencies in (15.31 15.5)) . If the regular g-polygon is constructible, 
the trigonometric functions in ( 15. 6 j) are algebraic numbers. This is the case for 

q E {2, 3, 4, 5, 6, 8, 10, 12, 15, 16, 17, 20, 24, 30, 32, 34, 40, 48, 51, 60, 64, ...} , (5.8) 

[4"8] Q Due to ( 15 .6p the constants ln(rfj) with di ^ 1, the divisors of q, and logarithms of further 
algebraic numbers will occur. Theses are all transcendental numbers 
Let us list the first few of these relations, see also 



IE + *l> 

lE + ^P 

JE + ^P 

lE + 4> 

lE + ^P 



lE + 4> 
lE + 4> 
lE + 1p 



lE + ll> [ -= 



= M2) 




(5.9) 






(5.10) 






(5.11) 


= ?ln(2) + - 
4 v ; 8 




(5.12) 


4 v ; 8 




(5.13) 


= ^[ln(2)-ln(V5-iy 


, l wc , , ^25 + 10^ 

+ - In 5 H 7r 

4 v ; 50 


(5.14) 



V5 
10 

v^s r 

10 



ln(2) - \n(V5 - 1) 



ln(2) - \n(V5 



- ln(5) + - 



40 



ln(2) - ln(\/5 



10 



iln(2) + iln(3) + ^V3 
iln(2) + iln(3)-^V3, etc. 



4 M5) 



4ln(5) 



10-2>/5 



40 



1 - 



'25 + 10^5 



50 



-7T 



5 
5 



7T 



(5.15) 



7T 



(5.16) 
(5.17) 

(5.18) 

(5.19) 



5.1.2 Alternating Single Sums at w = 1 

The alternating sums at w = 1 have the representation 



E 

fc=0 



-1) 



Ik + m 



1 

21 



m + I 
21 



with the reflection relation 
_1)* (_i)* 



E 

fc=0 



nk + I nk + (n — I) 



1 

2n 



1 2 2n 



. 1 Z 
^2 + 2^ 



ib[l 

S 2n 



(5.20) 



(5.2i; 



See [49] for special values of the trigonometric functions occurring in 
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For the sums at w = 1 one obtains 
.1)* 



E 

k=0 
oo 

E 

k=0 

oo 

E 

fc=0 
oo 

E 

fc=0 
oo 

E 

fe=0 
oo 

E 

fc=0 

oo 

E 

fe=0 
oo 

E 

fe=0 
oo 

E 

fe=0 
oo 

E 

fc=0 

oo 

E 

fc=0 



2k + 1 

(~l) fc 
3k + 1 

(~l) fc 
3k + 2 

(~l) fc 
4fc + 1 

(~l) fc 
4k + 3 

(~l) fc 
5fc + 1 

(~l) fe 
5fc + 2 

(~l) fc 
5fc + 3 

(~l) fc 
5A; + 4 

(~l) fc 
6A; + 1 

(~l) fe 
6fc + 5 



7T 

4 

1 

3 

1 

3 



7T 



+ ln(2) 



- ln(2) 



2v/2 
1 

271 
1 



7T 



7T 



-ln(V2-l) 



+ ln( v / 2- 1) 



7T 



1 + v^l ln(2) - -1= ln(V5 - 1) + * + = 

J V5 5^10 + 2^ 

1 - v^l ln(2) - 4= Ib(V5 - 1) + / V ^~ 1 

J V5 5V10 - 2v/5 

1 - v^l ln(2) + 4= ln(V5 - 1) + ~ 1 

J V5 5^10-2^ 

1 + V51 ln(2) + -1= ln(V5 - 1) + _ i + ^L 
J V5 5^10 + 2^ 



7T 



7T 



7T 



7T 1 

6 + ^ L 2 



- ln(2) - ln( y/3 - 1) 



7T 1 

6 



-ln(2)-ln(v / 3-l) 



etc. 



5.22) 
5.23) 
5.24) 
5.25) 
5.26) 
5.27) 
5.28) 
5.29) 
5.30) 
5.31) 
5.32) 



5.1.3 Single Infinite Sums of Higher Weight 

These sums obey representations which are obtained by repeated differentiation of (15. 3p and 
HQS} for m : 



E 

k=0 

00 

E 

fc=0 



1 1 1 ( m 

(ik + m) n F CH v n 'T 



(~l) fc 
(//c + m) n 



n-l 



T(n) (2/)™ 



r( 



(5.33) 
(5.34) 



Here Ch is the Hurwitz (^-function [52] with the serial representation 

00 



n=0 



(n + a) 4 



(5.35) 
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We consider the representations of the polygamma functions at rational arguments for p/q, 
p\q,qe N+, q < 120 One obtf ' 



ams : 
1 



(0 



with 



C: 



2n 



1)< +1 Z! (2' +1 - 1) + i , 



_1) 2(2^!^' 



where _B n denote the Bernoulli numbers [5U|55]. They are generated by 



x 



B 2n+1 = for n > 1 



exp(.x) — 1 



(5.36) 



(5.37) 



(5.38) 



n=0 



For odd values of Z no new basis elements occur due to (|5.36p . 
The reflection formula |36| 



^ (n) (! - z ) 



-ir 



ll) n {z) + 7T- COt(7Tz) 

v ' dz n y ' 



(5.39) 



implies the relations for the argument p/q and (q — p)/q. Likewise, the multiplication formula 
for the Hurwitz zeta function [52] holds, 



Ch(s, kz) 
rn l+1 ip^(mz) 



1 \ - . / n 

ks^u[s,Z + - 
n=0 

h— n \ / 



fc e N_ 



(5.40) 
(5.41) 



with the special relation 



m—l 



(-l)"' +1 Cn + l " 1] = 

fc=l 



(5.42) 



-,(0 



One obtains, cf. also 

^ (f) = (3^-l)Cw-^ 
For even values of Z, V (/) (l/3) and if)® (2/3) are linear in 7r' +1 -\/3 and C/+ij with 



--V3tt 3 -26C 3 

!V37r 3 -26C 3 
1 6 

V^tt 5 - 2904C 5 

o 

16 

— VStt 5 - 2904C 5 , etc. 
3 



(5.43) 

(5.44) 
(5.45) 
(5.46) 
(5.47) 



Special examples were also considered in 
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Therefore, only for odd values of I one new basis element due to ?/>( 2Z+1 )(l/3) contributes. 
The two values at argument 1/4 and 3/4 are related, see e.g. [57] . 



{21-1) 



\H-\ 



21 
4 a-i 



(2/) 



[^(2"- l)\B 2l \ +2(2/)!)8 D (20] 
— [7r 2 '(2 2 '-l)| J B 2Z |-2(2/)!/3 D (2/)] 
_ 2 a-i [ vr ^i| jE ; 2i | + 2(2/)!(2 2 ' +1 -l)C 2m ] 
+2 a-i [ 7r a+i|^|_2(20!(2 !B+1 -l)C a+ i] 



(5.48) 
(5.49) 
(5.50) 
(5.51) 



Here, .E^ denote the Euler numbers [58H61] . They are generated by 



E 



exp(:r) + exp(— x) 
f3 D is the Dirichlet /3-function [521162] 



/ — ' n! 



2n+l 



for n > 1 . 



(5.52) 



n=0 



°° C-T) fe 
fc=0 V y 



(5.53) 



which is also given by the inverse tangent integral Tii(x) [63,64j, being related to Clausen 
integrals [55] ■ The special value for I = 2 yields Catalan's constant 



~ J_l)fc 



C . 



(5.54) 



At even values of I no new constants appear. Odd values contribute with Ti(2Z). 
For z = 1/5,2/5,3/5,4/5 the relations 



/!(-l) m Cm [5' +1 - 1] 



(-1)' 
(-1)' 



^ (/) ' i) + ^ cotM L = 2 / 5 



I 2 = 1/5 



,(o ( i) + , w (n + ,<o (i U»g 



(5.55) 
(5.56) 
(5.57) 



hold. For even values of /, ip^(2/5) is thus dependent on ^>W(l/5). 
For 2 = 1/6, 5/6 the relation 



(I) 



[6 m - 3 m - 2 l+1 + 1] - V W 



(5.58 



holds. The reflection formula (I5.39P relates ^^-'(1/6) and ip^(5/6) for even values of I to linear 
combinations of 7r l+1 \/3 and (21+1- For odd I the -^-values can be expressed by ip^ (1/6)- The 
shift-formula by 1 /2 



2 l+1 if>®(2z) = il>®(z)+if)® (z+- 



(5.59) 
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implies 



if, 



(~) = (2 i+1 + 1) g) + (-I)'/! (3' +1 - 1) + i 



(5.60) 



and relates all values to ip® (1/3). 

For -2 = 1/8, 3/8, 5/8, 7/8 one obtains 



^ ( + ^ ( I) + 1>® { I) + 1>® ( I) = /K-l) m 0+i4 m (2' +1 - 1) 



Due to (15.411) and (I5.39j) the corresponding ■?/>- values can be expressed by ip® (1/8) 



1>® ( - 



(-1)' 



^ <o G) +7r ^ c ° tML=i / 8 

For z = 1/10, 3/10, 7/10, 9/10 two reflection relations and 
1>® Qq) + V> <0 (^) + 1>® (Jq) + V> W (^) = /!(-l)' +1 + i [10' +1 - 5< +1 - 2^ + 1] 



(5.61) 

(5.62) 
(5.63) 
(5.64) 



1 

To 

3 

To 

7 

To 

9 

To 



hold. By the shift relation ( I5.59P one obtains 



2'+yo(0_ v ,(o g 

2^V®(|)-^ g 



(5.65) 

(5.66) 
(5.67) 
(5.68) 
(5.69) 



No new constants contribute. 

For z — 1/12, 5/12, 7/12, 11/12 two reflection relations and 

1 \ m / 5 \ / 7 \ /ll 



12 



12 



12 



12 



«(-l)' +1 G +1 2 



x (6 /+1 - 3 m - 2 l+1 + 1) (5.70) 



hold. All values can be expressed by '0^(1/12) 
5 



12 
7 

T~2 



•1)' 



V 12 



2^0) (J) a 



^'5 



(-iy 



2 = 1/12 



etc. 



(5.71) 
(5.72) 
(5.73) 



27 



Eq. (Oljl yields 



12 



For odd values of I one obtains 



12 



lz=7/12 



(5.74) 



(5.75) 



Eq. (I5.70p does not imply a further relation. ^^(1/12) contributes as new constant for even 
values of I. 

For even values of I = 2q the Hurwitz ^-function (15.51) obeys the representation 



£ 

fc=0 



1 



(2qk + 2p- l) n (2q) n 



Ch n 



2p-l 
2g 



1 q 

-E 



k=l 



C n ( - I cos 



(2j9 - l)Jfe7T 



+ 5 n ( ^ ) sin 



(2p- l)jfe7T 



where CV and are represented by the Legendre unction [6 TJ (I5.77p . 



(5.76) 



(5.77) 



with@ 



C„(x) = Re;\> (exp(mx)) 
S^x) = Imx^ (exp(z7rx)) . 

For v G N, (I5.78P can be represented by Euler polynomials [HE] and powers of 7r, 

y 



5. 



C*2n ( - 
P 

2n+l 

' g 



4(2n- 1)! 

f-ir 



with 



4(2n)! 

2 exp(xt) 
exp(t) + 1 



7r 2 " +1 £;. 



P 

2n I - 



p,g G N+,p < g, 



(5.78) 

(5.79) 
(5.80) 



5> 

n=0 



n! 



(5.81) 



Since 



(m— 1) 



(l) = ( - 1)mr( — 1)! |pT 



(5.82) 



in particular xp^(l/8) and ^^(1/12), Z G N + ,/ > 1 consist of one term containing an integer 
power of 71 and a second term oc Re(\m)(xi+i(r)),r G Q according to (I5.76[ I5TTB1) . 



°We corrected typos in Eq. (10) of Ref. 
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In conclusion, the representation of the single sums of weight w > 2 and I < 6 require the 
additional constants : 

W, ^ (2fc+1) ( I) , Ti 2fc (l), tf« ( I) , ^ ( I) , ^ ( I) , ^ (±) , k G N + . (5.83) 



3 J V 5 / V 5 / V 8 / V 12 , 

Finally we would like to comment on a relation in Ramanujan's notebooks (69] Chapter 9, 
(11.3), which was claimed to involve a cyclotomic harmonic sum, 

oo r 

m = iS(MF-wlS(^TTF- (6 ' 85) 

with 

G(l) = . (5.86) 

This relation is unfortunately incorrect, cf. [69], p. 257. G(l) may be given by the following 
representation in multiple zeta values introducing nested harmonic sums [70] 

~ 0.16227193947148339072 (5.87) 

using the representations in [UJ. At present, no relation between the special constants used in 
(I5.87P are known. On the other hand, H(l) is given by 



mi) = x 



+ ^ i2) ( I 



2048 

0.14402290986880995023 . (5.88) 



5.2 Infinite multiple sums at higher depth 

As a further extension of the infinite cyclotomic harmonic sums [11] we consider the iterated 
summation of the terms (I4.60p . The corresponding sums diverge if the first indices have the 
pattern c« = 1, Sj = 1, ( 12. ip . However, these divergences can be regulated by polynomials in 
Co and cyclotomic harmonic sums, which are convergent for iV — > oo, very similar to the case 
of the usual harmonic sums [UIB]. We study the relations given in Section 14731 supplemented 
by those of the shuffle algebra (SH) of the cyclotomic harmonic polylogarithms (at argument 
x = 1), Eq. ( I3.2ip . In Table 3 we present the number of basis elements obtained applying the 
respective relations up to weight w = 6. The representation of all sums were computed by means 
of computer algebra in explicit form. We derive the cumulative basis, quoting only the new 
elements in the next weight. Up to w = 4 we derive also suitable integral representations over 
known functions. One possible choice of basis elements is : 

w = 1: 

^{i,o,i} = °o (5-89) 

^{1,0,-1} = - ln ( 2 ) ( 5 - 9 °) 

*{2,i,-i} = -1 + 7 ( 5 - 91 ) 
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weight 




A 


bn 


A cu 


A + sn + Mi 


A + bn + Mi + M 2 


A i C U i U i U i IV /I 

A + on + Mi + M2 + IVI 


1 


4 


A 

4 


4 


A 

4 


A 

4 


3 


3 


o 
Z 


on 
21) 


1 n 
10 


13 


3 


9 

3 


o 
z 


1 


Q 
O 


1 no 






ft 

u 


u 




Q 
O 


4 


500 


150 


163 


10 


10 


9 


6 


5 


2500 


624 


650 


21 


21 


19 


13 


6 


12500 


2580 


2635 


36 


36 


34 


25 



Table 3: Basis representations of the infinite cyclotomic harmonic sums over the alphabet 
{(±l) k /k, (±l) k /(2k + 1)} after applying the stuffle (A), shuffle (SH) relations, their combination, 
and their application together with the three multiple argument relations (Hi, H 2 , M), as far these 
lead not to new quantities. In the latter case we quote the cumulative number of basis elements 
appearing at the new weight. 



w = 2: 
w = 3: 



w = 4: 



<*{2,l,-2} 

^{1,0,3} 
C{1,0,-2},{2,1,-1} 

C{2,1,-2},{1,0,-1} 

C{1,0,-1},{1,0,1},{1,0,1},{1,0,1} 
^{2,1,-4} 
C{2,1 -3},{2,1 -1} 

C{1,0-2},{1,0-1},{2,1,-1} 



C{1,0-2},{2,1 -1},{1,0,1} 



-1 + C 



= Ca 



_______ I 

12 48 + 2 



= -Cln(2)+ I dx- 
Jo ' 



dx ^~ Li 2 (:r) 
o x + 1 

1 X2{Vx) 



X wx 



— Lu 



-1 -iXi{i) 

-(^M2)-| 6 ) (i-: 



1 



x 



2 J l + x 



1 



7T 



U 2 (-x) + — ) ln(l - x) 



+ ^Si i2 (a; 2 ) - Si j2 (x) - Si )2 (-x) 



(5.92) 

(5.93) 
(5.94) 

(5.95) 

(5.96) 
(5.97) 
(5.98) 



(5.99) 



~\ f dx-^- [In (l-y/x)- In (l + yfi)] Li 2 (x) 
z 7o l + x 

-ln(2) / 1 rf a ;-^Li 2 ( a ;) + ^[ln(2)7r-2C] (5.100) 
jg 1 + x z4 



0"{1,O-2},{2,1 -1},{2,1,1} - 



dx- 



x 



l + x 



ln(l — x)Li 2 (x) 



-i[l-ln(2)] fdx^-U 2 {x) 
Z Jo 1 -t- X 



30 



+ 



7T 



2 r 



24 



1 - - ln(2) 



C 



Here Si^x) denotes a Nielsen integral [7T] . 



- / — In 2 1 



(5.101) 



(5.102) 



At weight w = 5,6 we give only a few integral representations. They can in general be 
obtained form the Mellin transforms setting the kernel x N — > 1. The following basis elements 
are obtained : 



w = 5: 

°"{1,0,5} 



0"{1,0-1},{1,0,1},{1,0,1},{1,0,1},{1,0,1} 



Lis 



C{1,Q,-4},{2,1,-1} 



cr {l,0,4},{2,l,-l} 



7 
720' 



-It 



7 5 1 

2880 2 



dx- 



x 



1 + x 



Lh(x) 



7T 7T 

~90 + 360 



' 1 
+ - 



dx- 



LL 



-x) 



(5.103) 
(5.104) 

(5.105) 

(5.106) 



1 + x 

°"{2, 1,-4}, {1,0,-1}) °"{1,0,-3},{1,0,-1},{2,1,1}? °"{1,0,-3},{2,1,-1},{2,1,-1}? C{1,0,3},{2,1,-1},{2,1,-1} , 
°"{2, 1,-3}, {2, 1,-1}, {2, 1,1}? ^{l^, -2}, {1,0,-1}, {1,0,-1}, {2, 1,-1}? °"{1,0, -2}, {1,0,-1}, {2, 1,-1}, {1,0, -1} ? 
O"{l,0,-2}, {2,1,-1}, {1,0,1},{1,0,1}, 0"{1,0,-2},{2,1,-1},{2,1,1}, {1,0,1} • (5.107) 



w = 6: 

O"{l,0, -5}, {1,0,-1} = C-5,-l> (5.108) 

o"{i,o,-i},{i,o,i},{i,o,i},{i,o,i},{i,o,i}, {1,0,1} = -Li6 (^j , (5.109) 

o-{2,i,-6} = -1 - iXe{i), (5.110) 



^"{l.O, -2}, {2, 1,-1}, {2, 1,1}, {1,0,1}, {1,0,1}? ^{lA -2}, {2, 1,-1}, {1,0,1}, {1,0,1}, {1,0,1}? 
°"{1,0, -2}, {1,0,-1}, {1,0,-1}, {2, 1,-1}, {2, 1,1}? O"{l,0, -2}, {1,0,-1}, {1,0,-1}, {2, 1,-1}, {1,0,1}? 
a {l,0, -2}, {1,0,-1}, {1,0,-1}, {1,0,-1}, {2, 1,-1}? ^{lfl^}, {2, 1,-1}, {2, 1,-1}, {1,0,1}? 

°"{1,0, -3}, {1,0,-1}, {2, 1,1}, {1,0,1}? °"{l,0,-3}, {1,0,-1}, {1,0,1}, {2, 1,1}? cr {2, 1,-3}, {2, 1,-1}, {2, 1,-1}, {2, 1,-1} ? 
^{1,0, -3}, {2, 1,-1}, {2, 1,-1}, {1,0,-1}? °"{l,0,-3}, {2, 1,-1}, {1,0,-1}, {2, 1,-1}? O"{l,0,-3}, {1,0,-1}, {2, 1,-1}, {2, 1,-1} ? 
°"{2, 1,4}, {1,0,-1}, {2, 1,-1}? C{1,0,4}, {2, 1,-1}, {1,0,-1}? ^{l.O^}, {1,0,-1}, {2, 1, -1} ? 
°"{2, 1,-4}, {1,0,-1}, {2, 1,1}? 0"{2, 1,-4}, {1,0,-1}, {1,0,1}? °"{1,0,-4},{2,1,-1},{2,1,1}? 
^"{lA -4}, {2, 1,-1}, {1,0,1}? C{1,0, -4}, {1,0,-1}, {2, 1,-1}? °"{2,l,-4},{2,l,-2}? 

C"{2,1, -5}, {2,1,-1} • (5.111) 



Recently, infinite sums of a type proposed in 




were studied in [75], Eqs. (lla-c). 



They can be expressed in terms of sums studied in this Section : 

rrtS 1 bl>(k+l/2) —-0(1/2)1 /r110 , 
Jl{r) = 2L (2k + lY = ff {2,l,r} -^{2,l,r+l} + ^{2,l,r+l},{2,l,l} (5.112) 



7 For similar sums see [73] . 
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rr I \ 1 ^ Mk + 1/2) - V>(l/2)] 1 

^2(r) = -2^ 7^7- = ^{w},{2,-i,i} (5-113) 

fc=0 1 ' 

Moreover, values of usual harmonic polylogarithms [26] at x — 1, J c?t[Li p (±t) — Li p (±l)]/(l±t), 
are discussed. In part these integrals refer to all three letters of the corresponding alphabet. The 
corresponding representations involve infinite harmonic sums of depth d > 1 naturally, as e.g. 
for 

/ dJ^M = tf_ 1A0A1 (l) = -^C 5 ln(2) + -_ B| _ 1 (5.115) 
and in similar cases of higher weight, see Ref. [8|[TT||26]. 



5.3 Infinite multiple sums with more cyclotomic letters 

Let us now consider more cyclotomic letters. We study the sums up to weight w = 2 and 
cyclotomy I = 20 □, based on the sets of the non-alternating and alternating sums using the 
letters 



,,)/ w 1 < n < 2, 1 < I < 20, m < I . 
(Ik + m) n 



(5.116) 



We use the stuffle (quasi-shuffle) relations for the sums, the shuffle relations on the side of the 
associated cyclotomic harmonic polylogarithms, and the multiple argument relations for these 
sums, cf. Section 14.31 In some cases the latter request to include sums which are outside the 
above pattern. In this case the corresponding relations are not accounted for. At w = 1 the 
respective numbers of basis elements is summarized in Table 4. 



1 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


sums 


2 


4 


6 


8 


10 


12 


14 


16 


18 


20 


22 


24 


26 


28 


30 


32 


34 


36 


38 


40 


basis 


2 


3 


4 


5 


6 


6 


8 


9 


8 


10 


12 


10 


14 


14 


11 


17 


18 


14 


20 


18 


new basis 
sums 


2 


1 


2 


2 


4 


1 


6 


4 


4 


3 


10 


2 


12 


5 


3 


8 


16 


4 


18 


6 



Table 4: The number of the w = 1 cyclotomic harmonic sums (15.1160 up to I 
elements at fixed value of I, and the new basis elements in ascending sequence. 



20, the basis 



The reflection relation (I5.39P for the -^-functions for x ■<-»• (1 — x) implies that there are at most 
I + 1 basis elements. We showed that the use of the above analytic representations and the 
shuffle, stuffle, and multiple argument relations lead to the same number of basis elements for 
I < 6 in the non-alternating and alternating case. 
The independent sums at w = 1 up to I = 6 are : 

°"{1,0,1}' °"{1,0 ,-l}> °"{2,l,-l} ! °"{3,1,1}> C{3,1 ,-1}, <7{ 4) i CT{4,3-1}, C{5,1,1}) C{5,1,-1}> 

0"{5,2 -1}, 0"{5,3 -1}, C{6,1 -1} , (5.117) 

8 Relations between colored nested infinite harmonic sums have been investigated also in Refs. [76 |l77j recently. 



32 



see Sections 15.1.14 15.1.21 for equivalent representations. 
The dependent sums up to I = 6 are 



a {2,l,l} 
0"{3,2,1} 

°"{3,2,-l} 
a {4,l,l} 
a {4,3,l} 
a {5,2,l} 
°"{5,3,1} 
a {5,4,l} 

°"{5,4-l} 

°"{6,1,1} 
a {6,5,l} 
^{6,5,-1} 



1 - 0"{i,o,-i} + 2°"{iAl} 



O"{l,0,-1} — C{3,1,-1} + °"{3,1,1} 



1 1 

~2 ~ 3 

1 2 

2 + 2°"{i,o,-i} + cr{3,i ,-i} 



1 3 

2 ~ 4 
5 3 



- g - 4 CT {iA-i} + 4°"{W} - 2^2,1,-1} 
1 

^{1,0,-1} + ^{S.l.l} _ °"{5,2 -1} 

1 1 , 

— g - gCT{l,0,-l} - °"{5,l,-l) + °"{5,1,1} 

7 2 

12 ~ 
7 4 

r 

1 1 1 

-gO-{i,o,-i} + 2 a {i,i,-±} + ^^i 3 ' 1 ' 1 } 

7 2 1 

~10 ~ 3 <7 { 1 '°«- 1 -' ~~ a &> 1 -V + 2 CT{3 ' 1 ' 1} 
2 4 

+ 3°"{2,i -1} - °"{6,i -l}, etc. 



— - g^l.O-l} ~~ °"{5,1 -1} + cr {5,l,l} ~~ °"{S,3,-1} 



! g + gCT{l,0,-l} + °"{5,1,-1} - 0"{5,2,-l} + CT{5,3 -1} 



(5.118) 

(5.119) 
(5.120) 
(5.121) 
(5.122) 
(5.123) 
(5.124) 
(5.125) 
(5.126) 
(5.127) 
(5.128) 
(5.129) 



The remaining sums are related to those given in (I5.117l45.129p . The following counting relations 
for the basis elements were tested up to I = 700 using computer algebra methods. Let p,Pi,q be 
pairwise distinct primes > 2, and let k, hi be positive integers. The number of basis elements at 
w = 1 and cyclotomy I are given by 



r i+h 

( P - ly- 1 + 2, 

2y? (2 fc - i nr=ip9-n-i, 

(nr=i^)-^-2)-g+3, 

i ^(g fe - i nr=i^)-(^+2)(?-i), 



1 or I = 2 k 

k 



P 



2 fe n^i^ 



(5.130) 



Let us now consider the case w = 2. Applying the relations given in Section 14.31 and the 
shuffle algebra of the cyclotomic harmonic polylogarithms at argument x — 1 the results given 
in Table 5 are obtained for the number of basis elements. Again we solved the corresponding 
linear systems using computer algebra methods and derived the representations for the dependent 
sums analytically. 

The number of the weight w = 2 infinite sums for cyclotomy I is 



N s = 2/(2/ + 1) . 



(5.131) 



One may guess, based on the results for I < 20, counting relations for the length of the bases 
listed in Table 5. We found for all but the last column : 



N A (/) 



1(21 + 1) 



(5.132) 
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1 


Nc 

1 y b 


SH 


A 


A + SH 


A + SH + Hi 


A + SH + Hi + H 2 


A + SH + H x + H 2 + M 


1 


6 


4 


3 


1 


1 


1 


1 


2 


20 


13 


10 


3 


3 


2 


1 


3 


42 


27 


21 


7 


6 


6 


5 


4 


72 


46 


36 


12 


11 


10 


3 


5 


110 


70 


55 


19 


17 


17 


16 


6 


156 


99 


78 


27 


25 


24 


5 


7 


210 


133 


105 


37 


34 


34 


33 


8 


272 


172 


136 


48 


45 


44 


12 


9 


342 


216 


171 


61 


57 


57 


52 


10 


420 


265 


210 


75 


71 


70 


22 


11 


506 


319 


253 


91 


86 


86 


85 


12 


600 


378 


300 


108 


103 


102 


21 


13 


702 


442 


351 


127 


121 


121 


120 


14 


812 


551 


406 


147 


141 


140 


49 


15 


930 


585 


465 


169 


162 


162 


145 


16 


1056 


664 


528 


192 


185 


184 


50 


17 


1190 


748 


595 


217 


209 


209 


208 


18 


1332 


837 


666 


243 


235 


234 


63 


19 


1482 


931 


741 


271 


262 


262 


261 


20 


1640 


1030 


820 


300 


291 


290 


74 



Table 5: Number of basis elements of the w = 2 cyclotomic harmonic sums (15.1160 up to cyclotomy 
I = 20 after applying the quasi-shuffle algebra of the sums (A), the shuffle algebra of the cyclotomic 
harmonic polylogarithms (SH), and the three multiple argument relations (Hi, H2, M) of the sums. 



Nsh(0 
N a ,sh(0 
Na,sh,Hx(0 

Na,SH,Hi,H 2 (0 



(5/ + 3)/ 



6f 



l-(-l) 



6/ 2 



4/ + 7- (-r 



6/ 2 — 4Z + 3(1 



The latter relation fl 5 . 1 3 6 j) has been derived prior to us by D. Broadhurst ^ : 



Na.sh.Hx.h^O = + if(modp(l,2) = 0, 1, 3/4) 



and the corresponding generating function 

„3 



/(*) 



1 + 



X 



1 + x 



(5.133) 
(5.134) 
(5.135) 
(5.136) 

(5.137) 
(5.138) 



1=0 



We conjecture that in case of Na,sh,Hi,h 2 (0 the (M)-relations lead to a reduction of one in the 
basis for I being a prime. Otherwise quite significant reductions are obtained for which we do 
not know an explicit counting relation. The corresponding sequence is also not recorded yet in 
the data base 



5 We would like to thank D. Broadhurst for communicating this relation to us. 
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6 Generalized Harmonic Sums at Roots of Unity 

In Section [5] we considered real representations for the infinite cyclotomic harmonic sums. These 
are related to the infinite generalized harmonic sums at the roots of unity. We define 

lim S kl km (xi, ...,x m ; N) = cr kl km (xx, ...,x m ), (6.1) 

with S ku ___ ikm (xi, ■ ■■,x m ; N) a generalized harmonic sum (11.31) . see also [MHI], and Xj G C n , n > 1, 
with C n G {e„|e™ = 1, e n G C}; k\^ \ for X\ = 1. 

We seek the relations between the sums of w = 1,2. They can be expressed in terms of 
polylogarithms by : 

a w {x) = Li w (x), w G N, w > 1 (6.2) 
(Ji(x) = Lii(x) = — ln(l — x) (6.3) 

a 1A (x,y) = U 2 (x) + hn\l-x) + U 2 (-^^J (6.4) 

a lt i(x,x*) = Li 2 (a;) + iln 2 (l-a;) (6.5) 
Li w (a;) = U* w (x*) (6.6) 

and * denotes complex conjugation. Furthermore, the symmetric combination <j\ 7 \(x,y) + 
a hl (y,x) is given by 



0i,i (x, y) + 01,1(1/, x) = ln(l - x) ln(l -y) + U 2 (xy) . (6.7) 

Knowing the representations at w = 1 and the dilogarithms at the corresponding roots of unity 
one term <Ji t i(y, x) may be expressed by ( |6.7|) . 

In analyzing the functions 0fc 1 ,...,fc m (xi, ...,x m ) with x l k = 1 one may use the real representa- 
tions of nested cyclotomic sums for N — > 00, accounting for the respective sub-cycles at d, d\l 
and complex conjugation. 

An example is given by 

{ e i2 1 fcLi } = { e i2,e 6 ,e 4 ,e 3 ,el 2 ,e 2 ,e 5 12 *,e* 3 ,el,e* 6 ,e* 12 ,e 1 } . (6.8) 

Here we labeled the elements occurring in sub-cycles accordingly. 
The polylogarithms Lii(ef) and Li 2 (ef) obey : 

Im^ef)] = (6.9) 

Re[Li 2 (ef)] = ^"P + 'V . (6.10) 

More generally, PSLQ [79] tests let conjecture that Im [Li n (ef)l = r„ jZjfc 7r n for n odd and 
Re [Li n (ef )] = r n j jk n n for n odd with r n! i <k G Q. 

Now we extend Proposition 2.3 of Ref. [19], where we consider generalized harmonic sums 
S kl ,...,k m {xi,-..,x m ]N) with N G N, h G N + ,Xi G C, \xi\ < 1. Let / G N + and 

y\ = Xl (6.11) 

then 

m 

S kl _ km (x u ...,x m] N) = l[i k l S kl ,..., km (yi,..,y m ;lN) . (6.12) 
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Here the sum is over the Ith roots of Xj for i e Eq. (16.121) is called Distribution Relation. 

It follows from the Vieta's theorem [80] for ( 16. lip and properties of symmetric polynomials [81] . 
Eq. ( I6.12p contains the well-known duplication relation, cf. Eq. (2.15) [llj. 
If also x\ I for k\ = 1 the limit 

(Tku~,k m (xi, -,x m ) = lim S kl r .. }km (xi, ...,x m ;N) (6.13) 

N-¥oo 

exists. One may apply ( 16 . 1 2|I6 . 13j) to roots of unity X{ and yj, i.e. Xi = exp(2nini/mi) and 
Ujk = exp(2irikni/(mil)), k = I. ..(I — l),nj,mj G N+. Let us now consider the cases w = 1, 2 in 
more detail. 



6.1 w = 1 

The first element is real and occurs at I = 2 

Lii(e 2 ) = -ln(2), (6.14) 

representing the simplest alternating multiple zeta value, cf. e.g. [llj. At I = 3 we get the 
complex conjugate numbers 

I 77-7 

Li!(e 3 ) = --]n(3) + ^ (6.15) 

Lhiel) = -2 ln (3)--- (6.16) 

Due to (I6.9p . ln(3) and in are considered as basis elements from this level on. For all higher 
values of I one thus needs only to consider the real part of the w = 1 sums and one may work 
with the real representations given in Section [5j Let us consider the example 

Lii(e4) = — ln(l — e^j 



k=l k=l k=l k=l 

£ \4k-3 ~ 4k) + 62 ^ \4k-2 ~4k) +64 ^ \4k-l ~~ 4k) 

k=l v 7 fc=l v 7 fc=l v 7 



(6.17) 



Eq. ( I6.17P follows from 



n-l 



fc=i 

The type of sums occurring in (16.171) leads to digamma-functions and one may use their relations 
given before to find the corresponding basis representations. The first terms are given by : 

Li 1 (e 4 ) = --ln(2) + T (6.19) 
Lh(el) = Ii;(e 4 ) (6.20) 

Lii(e 5 ) = ^P^)-iM5) + i^ (6.21) 
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Lii(eD 

Re(Li 1 (e 5 )) + Re(Li 1 ( e 2)) 
Lii(e 6 ) 
Re(Li x (e 8 )) 
Re(Lii(ei2)) 



2 I 2 



- 2 M5) 

71) 

y 



i ln(5)+i-7T 



-iln(2)-iln(V2-l) 
- ln(2) - ln( V3 - 1) . 



In Table 6 we summarize the number of basis elements. 



(6.22) 

(6.23) 
(6.24) 
(6.25) 
(6.26) 
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9 
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11 
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2 


2 


3 


3 


4 


3 


4 


4 


6 


4 


7 


5 


6 


5 


9 


5 


10 


6 


Ref. 





1 


2 


2 


3 


3 


4 


3 


5 


4 


6 




















new 










































elements 





1 


2 





2 





3 


1 


2 





5 


1 


6 





2 


2 


8 





9 


2 



Table 6: The number of the basis elements spanning the w=l cyclotomic harmonic polylogarithms 
at Zth roots of unity up to 20. 



At cyclotomy I = 9 we find one basis element less than reported in [82]. The new elements 
contributing at the respective level of cyclotomy for I < 20 are : 



= 2 


ln(2) 


(6.27) 


= 3 


ln(3),vr 


(6.28) 


= 4 




(6.29) 


= 5 


Re(Li 1 (e 5 )),Re(Li 1 ( e 2)) 


(6.30) 


= 6 




(6.31) 


= 7 


Re(Li!(e?))| 3 fc=1 


(6.32) 


= 8 


Re(Li!(e 8 )) 


(6.33) 


= 9 


Re(Li 1 (e 9 )),Re(Li 1 ( e 2)) 


(6.34) 


= 10 




(6.35) 


= 11 


RetLi^))!^ 


(6.36) 


= 12 


Re(Li!(ei 2 )) 


(6.37) 


= 13 


ReCLiiCeaMjU 


(6.38) 


= 14 




(6.39) 


= 15 


Re(Li 1 (e 15 )),Re(Li 1 (e? 5 )) 


(6.40) 


= 16 


Re(Li 1 (e 16 )),Re(Li 1 (e? 6 )) 


(6.41) 


= 17 




(6.42) 


= 18 




(6.43) 


= 19 


Re(Li 1 (e 19 ))| 9 fe=1 


(6.44) 


= 20 


Re(Li 1 (e 20 )),Re(Li 1 (e| )) . 


(6.45) 
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6.2 w = 2 



We first consider the relations for Li 2 (e^). The following well-known representations for the 
function holds, cf. 183.1811 : 



Li 2 (e* d )=n 2 B 2 (J-\ +i C \ 2 (9) 



(6.46) 



with 



Ch(9) 
B 2 (x) 



E 

k=l 



sin(fc#) 



k 2 



1 oo 
7T 2 ^ 



k=l 



cos(27r/ca;) 



(6.47) 
(6.48) 



B 2 denotes the second modified Bernoulli polynomial. Due to (I6.10p only the imaginary parts 
have to be considered. The number ir occurs at w = 1, I = 3 only. The first terms are given by : 



Li 2 (ei) 
Li 2 (e 2 ) 

lm(Li 2 (e 3 )) 

lm(Li 2 (e 4 )) 

Li 2 (e^) 

lm(Li 2 (e 5 )) 



7T" 

if 

12 

Vs 

9 

c 

Li 2 (e 4 ) 



-7T 




7T 2 I 1 



7T 2 I 1 



V5 
1 

71 



(6.49) 
(6.50) 

(6.51) 

(6.52) 
(6.53) 



(6.54) 



lm(Li 2 (e 2 )) 



Li 2 (ej?) 
Li 2 (e^) 

lm(Li 2 (e 6 )) 

Li 2 (e^) 

lmLi 2 (e 8 ) 



lmLi 2 (ei 2 ) 



5V10 



1 

7E 



hG)-^ 1 ^)] 



i 

7E 



^(D f _ ^ i _ i 



i 

75 



Li* 2 (el) 

Li* 2 (e 5 ) 
3 

-lm(Li 2 (e 3 )) 
Li 2 (e 6 ) 

32 r 

24 v \3 



(i) + i (l _ 2 v2) c -i( 1+ v2; 



7T~ 



2 r ^ 

I 7T 

3 36 



(6.55) 

(6.56) 
(6.57) 

(6.58) 

(6.59) 

(6.60) 
(6.61) 
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The new basis elements spanning the dilogarithms of the Zth roots of unity for I < 20 are : 



= 1,2 




(6.62) 


1 = 3 


lm(Li 2 (e 3 )) 


(6.63) 


1 = 4 


C 


(6.64) 


1 = 5 


lm(Li 2 (e 5 )),lm(Li 2 (e 2 )) 


(6.65) 


1 = 6 




(6.66) 


1 = 7 


lm(Li 2 (e^))|' 


(6.67) 


1 = 8 


lm(Li 2 (e 8 )) 


(6.68) 


1 = 9 


lm(Li 2 (e 9 )),lm(Li 2 (e 2 )) 


(6.69) 


= 10 




(6.70) 


= 11 


lm(Li 2 (e*))l? , 


(6.71) 


= 12 


ImfLiofeio)) 


(6.72) 


= 13 


lm(Li 2 (ei 3 ))|? , 

V ^ V loy y \k=l 


(6.73) 


= 14 




(6.74) 


= 15 


ImfLiofeit;)) 


(6.75) 


= 16 


ImfLioCeiR^ ImfLiofeLH 


C6 76) 


= 17 


lm(Li 2 (ei 7 ))|fc = i 


(6.77) 


= 18 




(6.78) 


= 19 


lm(Li 2 (ei 9 ))|fc = i 


(6.79) 


= 20 


lm(Li 2 (e 20 )) . 


(6.80) 



Let us now turn to all convergent sums (J\^{x , y), x,y G C n . These sums belong to cyclotomy 
I if x = ef 1 , y = ef 2 and h,k 2 G N+, h < l,k 2 < I. 
We consider first the case x ^ 1, y — 1, 

(Ti,! (a;, 1) = Li 2 (x) + ^Li 2 (:r) 

= ^Re(Li2(x)) + ('r 2 -ir 1 2 ^7r 2 + 2[r 1 7rRe(Li 1 ( a ;)) + lm(Li 2 ( a ;))] , (6.81) 



with 



r\ = lm(— ln(l — x)) 
r 2 = Re(Li 2 (a;)) . 



(6.82) 
(6.83) 



<ri,i(e 2 ,x) = -i7r 2 + iln 2 (2) + Li 2 (i^ 



Including the basis elements of w = 1 up to I, no new basis element is obtained. Furthermore, 

(6.84) 

<ri,i(x,e 2 ) = -Li 1 (x)ln(2) + iLi 2 (x 2 )-Li 2 ( a :) + i [7r 2 + ln 2 (2)] - Li 2 {^^j ( 6 - 85 ) 

hold. Li 2 ((1 + x)/2) may yield a new basis element. In some cases besides x also — x is element 
of the the cycle of the roots of unity which have to be considered. Here, however, Li 2 ((1 + x)/2) 
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is given by 

Li 2 0^) = -Li 2 0^) + y - Lii(^)Lii(-^) - ln 2 (2) - ln(2) [Li x (x) + Li^-x)] . 



(6.86) 



Also the elements x = ejj and e± — x = 1 — x occur in the cycles, for which 



r 2 



Li 2 (l - e k n ) = -Li 2 (e*) + 2vrz-Li 1 (e^) + ^ (6.87) 

n 6 

holds. 

For I = 2 one obtains 

a ljl (e 2 ,l) = -ivr 2 + iln 2 (2) . (6.88) 

Because of ( I6.50P two basis elements contribute. If a corresponding special number has occured 
already at w = 1 at the same value of I, it is not counted as new. ir occurs at w = 2, I = 1, 
unlike for w = 1 where it contributes first at I = 3. No new basis element occurs at 1=2. For w 
= 3,4 the new basis elements occur for Li 2 (e^) only. We apply the relations, cf. e.g. [53] , 

Li 2 (— !— ) = Li 2 (-x)+ln(x)ln(l + x)-^ln 2 (l + x)+C 2 (6.89) 
\i + x J 2 

Li2 (l^) = -Li 2 Hr)-^ln 2 (l + o;) . (6.90) 
At I = 5 the above relations lead to corresponding reductions and the two functions 

Li2 , L , „ 91) 



\ ^5 / 

remain. For the first function the identity 



l + e 5 



Lb ( -^y? 1 ) = ^ (jf%) = - U *(-<) - \ ln2 (! + O (6-92) 
holds, through which the corresponding sum (T\^{e\ : e 5 ) can be expressed by 

es) = \^l{el) ~ U* 2 (e 5 ) - i (Li 2 (e 5 ))* + Li 1 (e 5 )*Li 1 (e 2 )* . (6.93) 



The second dilogarithm can be transformed in the following way : 



Furthermore, 



U 2 ( l \ e l + e f )=U 2 {-el) (6.95) 



holds, through which 



Li^-^-J = -^U 2 (e 5 ) + U;(el) + —n 2 --[U l (et)-U 1 (e 5 )} (6.96) 



40 



is obtained. 

The representations at w = 6 were given in [20], with Li4(l/2) the new basis element. The 
Clausen function Ci2(7r/3) used there is given by 

f \ [Li n (exp(-za;)) - Li n (exp(z'x))] , n even g 
1 | [Li n (exp(— ix)) + Li n (exp(zx))] , n odd 

Lsf (|) = §lm(Li 2 (e 3 )) (6.98) 

with 

Lsf\e) = -J dtt fc ln J '- fc - 

cf. [85]. In Table 7 we summarize the number of basis elements found at w = 2 using the relations 
in Section 4.3, the distribution relations (16.121 16.13p . and the relations for dilogarithms given 
above. We also list the number of basis elements for the class of dilogarithms at roots of unity, 
and in both cases the number of new elements beyond those being obtained at w = 1 at the 
same value of I. 
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8 





9 
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3 
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23 


40 


29 


Ref. [82] 


1 


1 


1 


1 


2 


2 


4 


4 


7 


6 


10 




















new 










































elements 
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1 


1 


2 


1 


4 


3 


5 


4 


10 


5 


14 


8 


12 


12 


24 


11 


30 


16 



Table 7: The number of the basis elements spanning the dilogarithms resp. w = 2 cyclotomic 
harmonic sums at /th roots of unity up to 20. 



Cl n (x) = 




(6.99) 



At w = 2 the respective new basis elements are : 



= 1 


7T 




(6.100) 


= 2 






(6.101) 


= 3 


lm(Li 2 (e 3 )) 




(6.102) 


= 4 


C 




(6.103) 


= 5 


lm(Li 2 (e 6 )),lm(Li 2 (ef)) 




(6.104) 


= 6 






(6.105) 


= 7 


lm(Li 2 (e£)) \l =l ,(Ti,i(e 7 ,e^) 




(6.106) 


= 8 


lm(Li 2 (e 8 )), a hl (e 8 , e 4 ), a ltl (e 8 , e 3 8 ) 




(6.107) 


= 9 


lm(Li 2 (e 9 )), lm(Li 2 (eg)), cri,i(e 9 , eg) 


,cJ"i,i( e 9,e 3 ), 0-1,1(69, e 3 ) 


(6.108) 


= 10 


ci,i(e 5 , e 2 ), cri,i{el, e 2 ), ai^eio, e 5 ) 


ci,i(eio, ei ) 


(6.109) 


= 11 


lm ( Li 2(en)) Li' <7i,i( e n> e ii) \ k=2 


°%i( e ii, e ii)L =3 , etc - 


(6.110) 



We mention that counting relations for majorants of the motivic numbers, which are claimed to 
be related to the bases of the sums crki,...,k m ( x ii •••) x m), %j £ C n , were given in [77]. The dimension 
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of the respective Q- vector space is majorized by the expansion coefficients D n (w) of 



oo 



5> n (w) 



X 



,n 



(6.111) 



n=0 



with 



1 



G 2 (a;) 



1 



1 — x — x 2 



G k (x) 



1 



k > 3 . 



(6.112) 



1- l[^(A;) + z/]x-(z/-l)x 2 ' 



Here, <f(k) denotes Euler's totient function [23] and v is the number of prime factors of k. 

7 Conclusions 

In evaluating massive 3-loop integrals with local operator insertions nested sums occur, contain- 
ing denominators which belong to the class of harmonic sums generated by cyclotomic polyno- 
mials. To simplify the computations, the relations between these quantities have to be know 
and used in computer algebra codes such as Sigma [17] and HarmonicSums [27J to allow for an 
efficient reduction in the corresponding summation problem. In integrals of similar structure we 
expect even more general terms (11.51) to occur. 

The usual harmonic sums [7J|8] and poly logarithms [26] were thus generalized to cover the 
newly occurring structures. We started with the harmonic polylogarithms extending the usual 
alphabet of denominators ranging to $2 0*0 1° general cyclotomic polynomials <& n (x), n > 3. The 
cyclotomic harmonic polylogarithms form a shuffle algebra. They have support x G [0, 1] and 
one may define a Mellin transform, usually of argument kN, with k,N G N + which span the 
finite cyclotomic harmonic sums (12.11) together with special values as the cyclotomic harmonic 
sums in the limit N — >■ 00. 

The cyclotomic harmonic sums are meromorphic functions with poles at the non-positive 
integers. They obey recurrence relations in terms of sums of lower depth and one may derive 
analytical asymptotic representations. In this way the cyclotomic harmonic sums are analytically 
continued from integer values of the sum index iV to iV G C. The cyclotomic harmonic sums form 
a quasi-shuffle algebra, cf . [9] . After the analytic continuation they obey differentiation relations, 
accounting for their values at iV — > 00. Furthermore, three multiple argument relations apply. 
Using these relations one may reduce the number of cyclotomic harmonic sums vastly growing 
with the weight to smaller bases. We study the case of the sums following from iteration of the 
denominators (I4.60p up to weight w = 5. Corresponding counting relations for the number of 
basis elements are obtained. 

The values of the cyclotomic harmonic sums for N — > 00, resp. the values of the cyclotomic 
harmonic polylogarithms at x — 1, are of special interest. They are linearly related to the infinite 
nested harmonic sums with roots of unity as numerator weight factors. Already in case of the 
single infinite cyclotomic harmonic sums a large number of new constants contribute beyond 
those know in the case of multiple zeta values and Euler-Zagier values [11] . The quasi-shuffle 
and multiple argument relations of the cyclotomic harmonic sums and the shuffle relations of 
the cyclotomic harmonic polylogarithms allow to derive basis representations induced by these 
relations. We studied in this respect the case of the infinite cyclotomic harmonic sums based 
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on the iteration of the summands (I4.60P to weight w = 6 and the sums of weight w = 1,2 for 
cyclotomy I < 20. Using computer algebra methods the explicit representation of all infinite 
cyclotomic harmonic sums were derived as well. For wide classes of relations explicit counting 
relations of the basis elements were given. The corresponding representations were derived with 
the Mathematica-based computer algebra system HarmonicSums [27J. 

The present investigations can be readily extended to cyclotomic harmonic sums and poly- 
logarithms of higher weight and cyclotomy, both in the case of finite values of N and for N — > oo, 
using the present methods. The requested computational time and storage resources grow ac- 
cordingly. This applies in particular to the derivation of the explicit representations of all sums 
over the corresponding bases. 

We finally considered also the generalized harmonic sums of weight w = 1,2 at Ith roots of 
unity for 1 < I < 20. They obey dilogarithmic reprentations. Besides the shuffle and stuffle 
relations, they obey distribution relations and the known relations for dilogarithms. We used 
these relations to derive corresponding basis representations. Compared to the case of the 
infinite cyclotomic harmonic sums these sums obey more symmetries. Thus at a given weight 
and cyclotomy I they are represented by a lower number of basis elements. We compared to 
results in literature. 
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A Appendix 

In the following we summarize some technical aspects needed to represent 
the present paper. 

First we summarize some aspects on cyclotomic polynomials, [23]. We 
tions of the polynomials 

x l + 1, le N\{0} 
in terms of cyclotomic polynomials for Z < 20 : 



For odd values of n, 



holds. The decomposition 



results thus into 



From 



it follows 



i i 

X T -L 






T 5 -1- 1 




Sf ^XJSf \0\X) 


x -pi 




sf^xjsfi^xj 


a; 7 + 1 


— 


$ 2 (x)$ l4 (x) 


r 9 4- 1 




Sf 2 l^X ) Sf 6 ^J^IS^J 


x 10 + l 


— 


$ 4 (x)$2 (x) 


x 11 + 1 




$ 2 (x)$ 22 (x) 


x 12 + 1 






x 13 + l 


= 


$ 2 (x)$ 26 (x) 


x 14 + l 




$ 4 (x)$ 28 (x) 


x 15 + l 




$ 2 (x)$ 6 (x)$io(x)$3o(a; 


x 17 + l 




$ 2 (x)$ 34 (x) 


.x 18 + 1 


— 


$ 4 (x)$ l2 (x)$ 36 (x) 


x 19 + l 




$ 2 (x)$ 39 (x) 


r 20 4- 1 










(x) = $„(-x) 


x 2k+1 


- 1 


= (x-l)JJ^(x) 

i 


x 2k+l 


+ 1 


= $ 2 (x) JJ$ 2i (x) . 

i 


| (x 2n - 


1) 


and $ 2 „(x) /(x™ 



expressions used in 
give the decomposi- 

(A.l) 



(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 
(A.IO) 
(A.ll) 
(A.12) 
(A.13) 
(A.14) 
(A.15) 
(A.16) 



(A.17) 



(A.18) 



(A.19) 



(A.20) 



§2n{x) | {X n + l) . 



(A.21) 
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If p is a prime and p\n then 



$Jx p ) 



(A.22) 



For n = 2 k , k G N + it follows that all $ 2 fc ( :E ) are cyclotomic polynomials. In (IA.6|A.12|A.14p 
more factors than just $2n(%) occur. They originate due to power-rescaling, i.e., 



x 15 + 1 



X 



15 



X 5 + 1 



y 5 + l 



y + 1 



My) ■ 



(A.23) 
(A.24) 



Therefore, all the factors of (x 5 + 1) and (x 3 + 1) have to emerge in the decomposition, and 
similarly for other N with more non-trivial divisors. For N = 2 k ■ n where the integer n > 1 is 
odd we get 



x 



2«.n + j 



X 2 " + 1 



y n + l 

y + i 



®2n{y) ■ 



(A.25) 



The argument remains valid if n is a product of odd primes. Therefore the only cyclotomic 
polynomials of the structure x a + 1 are those with a = 2 k , k G N+. 

For a proper definition of the cyclotomic harmonic polynomials which appear in sum repre- 
sentations like (12. 6p 



1 



x l ±r 



l G N_ 



(A.26) 



we perform partial fractioning. In the following we provide the corresponding expressions in 
terms of the words f l k (x) forming the cyclotomic harmonic polylogarithms up to / = 6 : 



(x 2 + l)- 1 



\ [/?(*) - fit?)] 

klfx 



X 



(X 



1) 



-1 



(x^-l)- 1 
(x 4 + l)- 1 
(x 5 - I)" 1 

(x 5 + I)" 1 
(x 6 - l)- 1 
(x 6 + I)" 1 



3 

f°s(x) 

1 



x 



2/ 1 (x)-/ 3 °(x)] 
/ 6 1 (x) + 2/ 6 °(x)] 
/ 2 °(x) - 2/ 4 °(x)] 



AV) - ~f&) - fax) - \fl{x) - \f!{x) 

/ 2 V) + \${X) - \fl{x) + ^/ 2 (X) - i/»(x) 



g [/i (*) - " 2/ s °(s) - tfC*) - 2/ 6 °(x) + 

i[/ 4 °(x) + 2A° 2 (x)-/ 2 2 (x)] etc. 



(A.27) 
(A.28) 

(A.29) 

(A.30) 

(A.31) 

(A.32) 

(A.33) 
(A.34) 

(A.35) 
(A.36) 

(A.37) 
(A.38) 
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B Appendix 



Proof of Eq. (14351, 



We proceed by induction on m. Let m — 1 : 

In j 2n 

^{ai.bi.cx} (2n) + 5' {ai)6l) _ ci} (2n) = ^ + ^ 



■1)* 



^(aii + fei) 01 ^(aii + 6i) ci 
y./ 1 1 



1 V(«i2« + &i) Cl (ai(2z - 1) + 6i) 
1 1 



+ 5J VK 2 ^ + &i) Cl (a 1 (2i-l) + 6 1 ) c i 

2 E (2ai , +6l)cl = 2 %^)W- 



In the following we use the abbreviation: 

M n ) := ^ ■S f {o m ,6 m ,±c m },...,{ai,6i,±ci> W- 

Now we assume that (I4.58P holds for m : 

^ ] >5'{a m+ i,6 m+ i,±c m+ i},...,{ai,6i,±ci} (27l) 



2n 



^ (a m +i« + m+ i) 6m + 1 ' 

2n (-lV 
+ S (a m+1 z + 6 m+1 )^ +1 S S K,W,*A,± C1 } (<) 

^ (a m +i« + & m +i) Cm+1 ^ (a m+ ii + & m +i) Cm+1 
A/ A(2i) A(2i - 1) 

I (On ... . ,i 4- /) , , Vm+i 



^ V (2am+i« + 6 m +i) Cm+1 ((2z - l)o m+ i + 6 m+ i) Cm +! 
l) 2i A(2i) (-l) 2i -M(2z - 1) 



i= i V (2a m+ ii + & m +i) Cm+1 ((2i - l)a m+ i + 6 m+ i) Cm + 1 
^ (2a m+1 i + 6™+!)°™+! 

n 1 

2 /_, 77^ ^~~~r v — 2 m S'{ 2a J c } {2oi,6i,ci} (*) 

f-j' (2a m+l i + m+ i) c '"+ i 

2 m+1 S'{2a m+1 ,6 m+ i,c m+ i},...,{2ai,6i,ci} (w) □ 
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Proof of Eq. (14331) . 



We proceed by induction on m. Let m — 1 : 

2n ^ 2n 

^{oLbi.ci} (2n) - S {ai , 6l _ ci} (2n) = ^ ^ + 



E 



^ (oil + &i) ci ^ (aii + 6i) Cl 
1 1 



^ V(ai2« + &i) Cl (a 1 (2i - 1) + hY 
1 1 



-E 



^ V(«i2« + &i) Cl (oi(2i-l)+6 



i) Cl 



n 1 
2 E pTlj^+^ = 2^2^-^} («) • 



In the following we use the abbreviation: 

A ( n ) '■= E dm " ' rf l 5 '{a™,6m,rfmC m },...,{a 1 ,fe 1 ,d 1 Cl} fa). 

Note that 



A(2n-1) = A(2n) - dm • ■ • di 



},...,{ai,6i,dici} (2n) 

(a m 2n + 6 m ) c ™ 

ic m _i},...,{ai,6i,(iici} 

(2n) 



(a m 2n + 6- 



A(2n) - ^d m _i • - -di 

. ^ ^ 7 j \ &m — 1 5 f*m — 1 1 — 1 Cm — 1 

},...,{ai,6i,dici} (2n) 

+ ^ ' ' ' dl + = • 4(2n) 



Now we assume that (14.591) holds for m : 

^ ] d m+ i ■ • • dl5{ am+li fr m+1 ,d m+ ic m+ i},...,{ai,&i,a!iei} (2fl) 



2n | 

/ . 7 ■ , t \ r ~ / , d m " " " diS{ a b m ,d m c m } ,—,{ai ,bi ,dm} (*) 

2n 



" E f„ a I, Vm+i E ^ m ' " ( ^l'S'{a m ,f) m ,d m c m },...,{ai,bi,dici} («) 



^ (a m+ ii + & m+ i) Cm+1 ^ (a m+ iz + & m +i) Cm+1 
^-v / A(2i) , A(2i-1) 



^ V (2a m+ ii + 6 m+ i) Cm +! ((2i - l)a m+ i + b m+1 ) c ^+ 1 
_ ys ( (-l) 2 Vl(2z) (-1) 2 -M(2* - 1) 



^ V (2a m+ ii + 6 m+ i) c ™+ 1 ((2z - l)a m+ i + b m+ i) c ™+^ 
2 y A(2i-1) A(2z) 



^ ((2z - l)a m+ i + 6 m+ i) c ™+! ^ ((2z - l)a m+1 + 6 m+ i) Cm+1 

n 1 

2 / 77^ r~ r ~ 2 5'{2a m ,6 m — a m ,c m },...,{2ai,6i— ai,ci} (^) 

^ (2a m+1 z + 6 m+ i - a m+1 ) c -+i 
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nm+l a f \ I — I 

z °{2a rn+1 ,b m+1 -a rn+1 ,c m+1 },...,{2a 1 ,b 1 ~a 1 ,c 1 } \ lb ) <— 1 



Proof of Eq. §£M\£M> 



We start with CfO^x). We consider integrals of the form 

rx i 00 j „ 2i+Cj 00 n i rx 



E(2i + c,)-+^» (<) 



Summing over j yields the desired result. 
For CfO^(x) we consider the integrals 



2^ 2i + c i + 3 2w 2 Jfc + c,.).* bW 
i=i j fc=i j 

= E (2»+ Cj . + i)^- M >- B(i) • 

Summing over j yields the desired result. The case Cfi^(x) follows analogously. 
For CfO^(x) we consider the integrals of the form 



(2i + c,Y 



~l —^w^r s " {i)iv 

rx 00 00 i+l„,2i+c 3 -+2 



48 



x oo i 2fc+c 3 -+2 

■>o (2fc + c,- + 2)- 

/•x 00 * _.fe+l ?/ 2fc+ Cj +2 

[ x f2y 2i+c * +1 S {2 , Cj ,„ a} , n (i + i)dy 

Jo i=0 

^ / ^^ +3 % C3 , CTfl} ,n(^ + l)^ 



i=0 

x 2i+c ? + l 



- XJ(2i + c . + i)^ 2 - c - CTa >'" W 

-IJ(2i + c . + 2)' 5{2 ' c - CTa} ' nW • 
Summing over j yields the desired result. The case CfO ^(x) follows analogously. 
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